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Dust Processing and Grain Growth in Protoplanetary Disks in 
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ABSTRACT 

Mid-infrared spectra of 65 T Tauri stars (TTS) taken with the Infrared Spec- 
trograph (IRS) on board the Spitzer Space Telescope are modeled using popula- 
tions of optically thin dust at two temperatures to probe the radial variation in 
dust composition in the uppermost layers of protoplanetary disks. Most spectra 
with narrow emission features associated with crystalline silicates require Mg- 
rich minerals and silica, but a very small number suggest other components. 
Spectra indicating large amounts of enstatite at higher temperatures (400-500 K) 
also require crystalline silicates, either enstatite or forsterite, at temperatures 
lower (100-200 K) than those required for spectra showing high abundance of 
other crystalline silicates. A few spectra show 10 fim complexes of very small 
equivalent width. They are fit well using abundant crystalline sihcates but very 
few large grains, inconsistent with the expectation that low peak-to-continuum 
ratio of the 10 fim complex always indicates grain growth. Most spectra in our 
sample are fit well without using the opacities of large crystalline silicate grains. 
If large grains grow by agglomeration of submicron grains of all dust types, the 
amorphous silicate components of these aggregates must typically be more abun- 
dant than the crystalline silicate components. We also find that the more there 
is of one crystalline dust species, the more there is of the others. This suggests 
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that crystalline silicates are processed directly from amorphous silicates, whether 
by evaporation of the amorphous grains and condensation in chemical equilib- 
rium or by annealing of the amorphous precursor grains. This also suggests that 
neither forsterite, enstatite, nor silica are intermediate steps along the way to 
producing either of the other two for the majority of the crystalline dust pro- 
duced. Crystalline silicate abundance is correlated tightly with disk geometry, 
in the sense of higher crystallinity accompanying more settled disks, which are 
commonly associated with growth and settling of grains. Large-grain abundance 
is also correlated with disks that are more highly settled, but with a wide range 
of large grain abundance for a given degree of settling. We interpret this range 
to mean that the settling of large grains is sensitive to individual disk properties. 
We also find that lower-mass stars have higher abundances of large grains in their 
inner regions. 

Subject headings: circumstellar matter, infrared: stars, stars: pre-main-sequence, 
planetary systems: protoplanetary disks 



Introduction 



Comets are thought to be largely unaltered reservoirs of leftover material from the pri- 



mord ial mixture of dust in the Solar Nebula orbiting the newly-formed S u n (e.g ., IWooden et al 



20051 ). Spectral observations of Comet Hale-Bopp 



frared Space Observatory {ISO; iKessler et al.l Il996) all o wed for detailed spectral models 



over a wide range of wavelengths (ILi &: Greenberg 



by 



Crovisier et al.l (119971 ) using the In 



1998; Wooden et al. 1999; Barker et al. 



20021 ) . The amorphous dust of pyroxene ([Mg,Fe]Si03) or olivine ([Mg,Fe]2Si04) stoichiom- 
etry (henceforth, "amorphous pyroxene" or "amorphous olivine", respectively) was found 
to require a stoichiometric ratio of magnesium to iron of about unity, being relatively 
i ron-rich (IHarker et al.l |2002| ): the crystalline pyroxene must be iron-poor, as it is cooler 
(IWooden et al.lll999l): and th e larger grain s mus t be porous, fluffy aggregates of smaller 



grains ( 



Bradley 



Li &: Greenbergl Il998l : iHarker et al.l |2002| ). Interplanetary Dust Particles (IDPs; 



2003h . believed to originate from comets, were often found to be highly porous 



aggregates of Mg-rich crystalline silicate dust, having qualities very similar to those inferred 
from comet spectra and spectra of Young Stellar Objects (YSOs; systems in the process of 
forming stars and planets). 

Using increasingly powerful telescopes and sensitive detectors, astr ophysicists have ob- 
tained increasingly higher-quality spectra for decreasingly fainter objects. ICohen fc Witteborn 
(119851 ) distinguished 10 fim silicate complexes in absorption versus emission, while iHonda et al. 
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( 120031 ) identified amorphous silicates, forsterite, enstatite, and silica using their spectrum of 
the 10 yum complex of Hen 3-600A. 

Approaches toward analysis of the dust composition of protoplanetary disks in Her- 
big Ae/Be and T Tauri (SED Class II YSO) systems have taken two forms: computa- 
tion of indices from flux ratios over the 10 fim (and sometimes longer wavelength) com- 
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T he Infrared Spectro graph (IRS; iHouck et al.l |2004| ) on board the Spitzer Space Tele- 
scope (IWerner et al.ll2004l ) allows good S/N spectra to be obtained at 5-37 /im wavelength for 



objects with flux densities of a few tens of mJy. Dust at a considerable range of temperatures 
in protoplanetary disks gives rise to the emission seen from s uch systems in IRS spectra. Our 



approach has been to mod el this using two temperatures (see lChen et al.ll2006l: iKastner et al. 



2006 : ISargent et al.ll2008l ). though others have used more temperatures (see iBouwman et al. 
20081 ) . Here we model the dust emission over the 7.7-37 /im w avelength range for 65 Class 
II YSOs in the 1-2 Myr old Taurus- Auriga star- forming region (IKenyon Sz Hartmannlll995l ) 
using dust at two temperatures, in order to analyze the dust composition in the inner and 
outer disk regions. 

"Crystalline indices" computed from flux ratios that characterize crystalline silicate 
emission have the advantage of being quickly computed and measured repeatably for a wide 
range of objects. However, they do not make full use of the information on dust emissivity 
available in a high-quality mid-infrared spectrum; for instance, they measure the combined 



effect of grain growth and crystallinity, but such are difficult to disentangle (jWatson et al. 
20081 ) without detailed modeling using dust opacities, as is done here. 



Data Reduction 



Much of the process of data r educt ion of the 65 s pectra analyzed in this study has been 
described by ISargent et al.l (120081 ) and iFurlan et al. (l2006h . Here we give a brief overview 
and provide details specific to this study. We note that the spectra of ZZ Tau and ROXs 42C 
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used here are exactly the same as used and described by lSargent et ahl (120081 ): the following 
discussion of data reduction applies to all spectra except those of ZZ Tau and ROXs 42C. 



2.1. Observations 

Sixty-three of the 65 ob jects in our s ample come from the 85 Class II Young Stellar 
Objects (YS Os) analyzed by iFurlan et al.l (120061 ) . The other 22 objects whose spectra were 
analyzed by Furlan et al. ( 20061 ) were rejected from our study due to high (here, defined as 
having Ay > 6) or unknown extinction, the presence of Polycyclic Aromatic Hydrocarbon 
bands, saturation of the detectors, an d confusion of Class II YSOs with Class I YSOs in the 
same spectrograph slit. To the 63 from lFurlan et al.l (120061 ) we add GG Tau B and Haro 6-28. 



For op acity calibration (see Section 3), we reduced data for ROXs 42C (see ISargent et al. 



( I2OO8I )). HBC 656, and TW Cha. All of these objects were observed using the combinations 
of IRS modules Short-Low and Long-Low or Short-Low, Sho rt-High, and Long-H igh in either 
Staring Mode or Mapping Mode; for more information, see lSargent et al.l (120081 ). 



2.2. Extraction and Calibration of Spectra and Rebinning of High Resolution 

Spectra 

We obtained Basic Calibrated Data products from the Spitzer Science Center (SSC) 
for each of the targets in our sample. BCD data are corrected for the dark current, stray 
light, and flatfield variations. The spectrum of HBC 656 was reduced using data from the 
S16.1.0 data calibration pipeline. The spectrum of ZZ Tau was reduce d using data from th e 
S14.0.0 IRS data calibration pipeline, using the methods described by lSargent et al.l (120081 ): 
the present discussion of extraction and calibration of spectra excludes the spectrum of ZZ 
Tau. For the spectra of TW Cha and the 64 T Tauri stars in Taurus-Auriga excluding ZZ 
Tau, the data are from the S13.2.0 data calibration pipeline. 

Next we identified arid fixe d by interpolation i n the spectral direction bad pixels, as 
described by lWatson et al.l (120071 ) and lSargent et al.l (120081 ). We th en extracted spectra using 
the Spectral Modeling, Analysis, and Reduction Tool (SMART; iHigdon et a3l2004[ ). We 
sky-subtracted data from the low spectral resolution modules SL and LL using the off-order 
observation. We did not subtra c t sky f rom high spectral resolution (SH and LH) observations, 
as explained by ISargent et al.l (120061 ). Multiple Data Collection Events (DCEs; detector 
samplings) at the same telecope position, same order, and same modules were averaged, and 
the uncertainties were calculated from the standard deviation of the mean. Tapered column 
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extractio n was used for SL a nd LL data, and full-column extraction was used for SH and 
LH data (ISargent et al.l 120061). We us ed Relative Spectral Response Functions (RSRFs) to 



calibrate the flux ( Sargent et al.l [20061 ). For both orders of SL and for the both orders of LL 
(for all observations except for HBC 656), a spectral template of a Lacertae (Al V; M. Cohen 
2004, private communication) of higher spectral resolution than the templates described by 
Cohen et al.l (120031 ) was used. For the LL flux calibration of HBC 656 over second, bonus, 
and first orders out to 36 /im wavelength, RSRFs generated from data of ^ Dra and the 
template for ^ Dra from ICohen et al.l (120031 ) were used. Past 36 /xm, a LL first-order RSRF 
generated from data of Markarian 231 and the template for M arkarian 231 (J. Marshall, 
private communication: iMarshall et al.ll2007: lArmus et al. 120071) were used. For SH and LH, 



the spectral template for ^ Dra (K2 HI) by ICohen et al.l (120031 ) and its IRS data were used 
to generate RSRFs. Except for GG Tau A, GG Tau B, ZZ Tau, RO Xs 42C, HBC 656, a nd 
TW Cha, all RSRF-calibrated spectra are the same as presented by lWatson et al.l (120081 ). 



We rebin our high-resolution spectra (see Table Al for whether the spectra were obtained 
with high or low spectral resolution) to low spectral resolution. The width of each bin 
is determined by requiring spectral resolution R=A/AA ~ 90, and rebinned spectrum is 
sampled such that each bin is AA from the next bin. 



2.3. 



Correction for Extinction, Mispointing, Uncertainties, and Stellar 

Photosphere 



We corrected our spectra for extinction using the same method described by lSargent et al 



( l2006l ). except we assume a ratio of V-band extinction to total opt i cal de pth at 9.7 /im, 
Ay/xg 7 of 25. This value is s u ggest ed by data plotted by lChiar et al.l (120071 ). We use values 
of Ay given by iFurlan et al.l (120061 ). We do not correct for extinctio n if Ay < 1.4, as w e 
want to minimize any artifacts due to overcorrection of extinction (see lSargent et al.ll2006l ). 
CoKu Tau/4, GG Tau A, GG Tau B, and Haro 6 -28 are not corrected, but V410 Anon 



13 is, assuming Ay = 5.8 (s ee ISargent et a. 



20061 ). The extinction correction for ROXs 



42 C is the same a s mad e by Sargent et al. ( 20081 ). The average (V-R) colors as measured 
bv iBatalha et all Jl998[ ) of HBC 656 (also known as AS 216) are ~ 0.64, which is slightly 



bluer than the intri nsic (V-R) color of a K2 star (the spectral type assigned to HBC 656 by 
Batalha et al.lll998l ) of 0.74 given by iKenyon fc HartmannI (119951 ). This is consistent with 
HBC 656 being of slightly earlier spectral type (e.g., KO) and having Ay ~ 0, so we assign 
Ay = to HB C 656. Finally, we com pute Ay = 1.78 for TW Cha using I and J magnitudes 
from DENIS (ICambresv et al.lll998l ) and intrinsic colors for a K4 main-sequence star from 
Kenyon fc HartmannI (119951 ). assigning a spectral type of K4 to TW Cha despite iLuhman 
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( 120041 ) adopting K8 for TW Cha after having found various classifications for this star in the 
hterature of KO, K8, and MO. 

Next, we corrected for shght mispointing of the telescope from the s tandard nod posi- 
tions by scaling the mispointed nods, as described by lSargent et al.l (120081 ) . The scalars used 
are given in Table Al and further described in the Appendix. After applying these scalars, 
the spectral uncertainties were computed. For all objects except ZZ Tau, the flux density 
unc ertainties were half t he difference between the spectra obtained at the two nod positions 



[ see iSargent et al.ll2006l ). The uncertainties were derived from only two independent mea- 
surements of flux (except for ZZ Tau), so sometimes the relative uncertainty was less than 
1%. Because of this and because of uncontrolled systema tic errors, when the r elative uncer- 
tainty was less than 1%, it was set to 1% as described by lSargent et al.l (120081 ). The spectra 
plotted by iFurlan et al.l (120061 ) for CoKu Tau/4, DM Tau, and GM Aur showed very little 
excess emission above the stellar photosphere at A < 8 fim; for these, we subtracted Planck 
functions using solid angles and effective stellar temperatures given by ISargent et al.l (120061 ) 
for the same purpose. P hotospheres w' e re sub tracted for ZZ Tau and ROXs 42C in exactly 
the manner described by lSargent et al.l (l2008l ). For all other spectra, the stellar photosphere 
contribution to the spectra appears minor, and stellar photospheres are not subtracted before 
modeling them. 



3. Analysis 
3.1. Models 

To characterize the silicate dust composition of the upper levels of the T T auri disks giv- 
ing r ise to the emission features in our spectra, we construct detailed models (ISargent et al. 
20081 ). Briefly summarized, each model is a sum of emission from blackbodies at tempera- 
tures T^ (warm) and T^ (cool) and emission from optically thin dust grains proportional to 
their opacities multiplied by Planck functions at T^„ and T^. 



(A) 

i 



+ 



where Bi,(A,T) is the Planck function evaluated at wavelength A and temperature T, Vt^ 
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is the solid angle of the blackbody at temperature Tc (T^), ac,i is the mass weight of 

component j at temperature Tc (T^), and Kj{X) is the opacity (effective cross-section per unit 
mass) at wavelength A of dust species j. The mass weights, a^j (a^^j), are equal to rric^i/cP 
[iTiw^i/cP), where m^j (jn^j) is the mass of dust species i at Tc {T^,) and d is the distance 
to the T Tauri star, which is assumed to be 1 40 parsecs for T Tauri stars in the Taurus- 
Auriga star- forming region (IKenyon et al.lll994l ). The blackbody components of each model, 
ficBi.(A,Tc) and f2^By(A,T^), represent optically thick continuum emission, unsubtracted 
stellar photosphere emission (Fjy proportional to 1/A^, and it should be small for spectra for 
which we have not already subtracted a Planck function at the stellar effective temperature 
to represent it), and emission from grains in the optically thin disk regions whose opacities 
are featureless continua, like amorphous carbon and very large silicate grains. The model 
components proportional to opacity multiplied by a Planck function represent emission from 
submicron- to few-micron-sized dust grains, located in the optically thin disk regions, with 
strong resonances giving rise to detectable features in their opacities. 

We will assume in the following that all grain species in the s ame location in the di s k have 
the same temperature, the thermal coupling provided by gas. iKamp fc DuUemondl (l2004l ) 
find that the layer with Ay < 0.1 (where Ay is the visual extinction of star into the disk), 
the layer above the layer where dust emission features arise (what they call the "superheated 
surface layer"), has gas and dust temperatures within 10% of each other, suggesting the gas 
would bring dust grains of differing optical propertie s to roughly the same temperature. In 
general, settling of dust towards the disk midplane (jP'Alessio et al.ll2006l ) should enhance 
the gas particle density for the optically thin regions in the disk atmosphere that gives rise to 
the 10 /xm silicate feature becau se the gas density is higher towards the disk midplane (see 
Figure 1 of lGlassgold et al.ll2004r). so dust and gas tempe ratures should be better coupled in 



more settled disks. However, IChiang fc GoldreichI (119971 ) found that, for a fiducial disk at 1 
AU, the dust grains should remain at their radiative equilibrium temperatures, undisturbed 
by gas-grain collisions, for the vast majority of plausible gas particle densities, 10* to 3x 10^^ 



14 



cm~^, at this distance from the star. They found that only between about 10^"^ and 3x 10 
cm~ ^ would the dust gra i n tera peratures be affected by collisions with gas molecules. At 8 
AU, iGorti fc HoUenbachl (120081 ) find that gas and dust temperatures are only equal for Ay 
slightly greater than 1, suggesting gr eater difference betw e en du st and gas temperature in 
the "superheated surface layer" than lKamp fc DuUemondl (12004 ) predict. For the purposes 
of dust modeling of T Tauri stars, Planck functions and opacity-weighted Planck functions 
at two temperatures (at most) will be used, while keeping in mind the results of these 
sophisticated studies of dust and gas coupling. We allow for the two temperatures in order 
to study compositional differences of disks between inner disk regions and outer disk regions. 
If, however, there is not sufficient couphng of dust and gas heating in the optically thin 
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surface layers giving rise to emission features in our TTS spectra, this means dust that 
absorbs starhght poorly at optical and near-infrared wavelengths with respect to mid-infrared 
wavelengths, like enstatite and forsterite, will be located closer to the star than dust like iron- 
rich amorphous pyroxene and olivine at the same temperature. Such iron-rich amorphous 
silicates absorb starlight more efficiently at visibl e and near-infrared wavelengths with respect 
to mid- infrared wavelengths (IHarker et al.ll2002l ). 

To determine the best fit of a model to a spectrum, we minimize between the model 
and the 7.7-37 /im region of the spectrum, with respect to the two temperatures, the solid 
angles, and the mass weights. We choose ranges of Tc and to explore, and, picking a 
single pair of [Tc,Ti„], we minimize x^; 



X 



(2) 



with respect to all solid angles and all mass weights. When this results in mass weights 
or solid angles that are negative, we eliminate from the model the component whose fiux 
integrated from 7.7-37 /im is most negative, and minimize of the model without this 
component. This is repeated at a given temperature pair [Tc,T^] until there are no negative 
mass weig htfl This entire process is carried out for each temperature pair. The temperature 
pair for which per degree of freedom (where the number of degrees of freedom is the 
number of data points (~ 200) minus the number of temperatures used, 2, minus the number 
of nonzero solid angles minus the number of nonzero mass weights) is lowest is the best-fit 
model for a spectrum. 

We also determine uncertainties for our model parameters. We define the uncertainty in 
a parameter to be the increment of the parameter from its best-fit model value that changes 
by 1.0 the value of per degree of freedom, or x^, of the model with the incremented value 
of the parameter. We use a Taylor series expansion to approximate this. 



We note that our definition of uncertainty in a parameter is similar to what lWooden et al 



(119991 ) call "independent point uncertainty" in that all model parameters are kept at their 
best-fit values except for the one whose uncertainty is being determined. However, our uncer- 
tainties are determined differently from how others have been determined in the literature. 
Press et al.l (Il992l ) suggest taking the uncertainty to be that increment of the parameter 



^The only solid angle that was found to be negative and in the subsequent iteration zeroed out was that 
of the cool blackbody component of CoKu Tau/3 



- 9 - 



that increases by a number c omputed from the incomplete gamma function for a given 
level of significance. lAvnil (119761 ) found that uncertainties computed using the increment in 
determined in this manner are quite consistent with uncertainties computed from Monte 
Carlo analysis. In such an analysis, the standard deviation in the set of parameter values 
obtained by fitting multiple versions of the same spectrum, each manipulated by adding 
noise according to Poisson statistics, is taken to be the uncertainty in that parameter. 

For example, the model for DG Tau has only eight constraints. The incomplete gamma 
function implies that an increment of of 9-31 would determine the 68.3% (one sigma) 
uncertainty level on parameters. This would result in a one sigma uncertainty on the mass 
weight of cool small amorphous pyroxene of 1.68 X 10"^*^ -2^. This is nearly identical to the 
uncertainty on this component obtained via the Monte Carlo analysis described previously 
of 1.62 X 10^^^ J^. Using the method adopted in this study, an uncertainty of 8.21 x 10^^^ 

is obtained, which is almost 5 times higher than the one sigma uncertainties obtained 

cm'' ' c3 o 

using the incomplete gamma function to determine the increment of x^ or the Monte Carlo 
approach. We note that, since the number of data points (about 200) in one of our spectra 
is much greater than the number of model constraints (by a factor of >20), o ur uncertainty 
level involves increasing x^ by much more than suggested by iPress et al.l (119921 ) . By so doing, 
we conservatively bias the uncertainties in our parameter values to many times one sigma. 



Our models include the opacity of annealed silica from iFabian et al.l (120001 ). This was 
shown to provide the best fit to emiss ion features at 9, 12, 16, and 20 /xm in astronomical 
spectra of TTS by lSargent et al.l (120081 ). We now justify our choice of opacities for submicron 
amorphous and crystalline silicates and large (few-micron-sized) silicate dust grains. In order 
to choose the best opacity curves to use to model emission from TTS, we assume that all 
disks in our sample have the same major kinds of dust, but in different relative abundances. 
Assuming all disks have basically the same major kinds of dust - amorphous silicates of olivine 
and pyroxene composition, crystalline olivines, crystalline pyroxenes, and crystalline silica 
- our goal is to find the best realistic opacities to represent these major kinds of dust, not 
necessarily to constrain the specific properties of the grains themselves (e.g., sizes, porosities, 
composition). 



3.2. Forsterite Opacity 



The first opacity for our models with which we concern ourselves is that of forsterite. 
Olivines have been used before in modeling dust emission to fit emis sion (and sonaetime s 
absorption) complexes at 10.0, 11.2, 16, 19, 23, 28, and 33 yum (e.g.. iMolster et al.l 120021 ). 
One spectrum for which these features appear at very high signal-to-noise ratio (S/N) is that 
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of HBC 656, which we include as our first and most important forsterite exemplar. Others 
with these features in their spectra at fairly high S/N ratios that we include as forsterite 
exemplars are DK Tau, GN Tau, IS Ta. u, ROXs 42C (which has already been identified as a 
silica exemplar; see lSargent et al.ll2008l ). V836 Tau, and V955 Tau. 



To fit the forsterite features in these exemplars' spectra, we have com piled opacities of 
crystalline olivine from a wide range of sources. From I Jaeger et al.l (119981 ). we include opac- 
ities from transmission measurements of the following olivine powders pressed into potas- 
sium bromide (KBr) pellets: forsterite (Mg2Si04), magnesium-rich olivine (Mgi gFco. 28104), 
hortonolite (MgFeSi04), and fayalite (Fe2Si04). To explore further the effects on opacity of 
the full range of iron content in the olivine solid series (olivines of all Mg/Fe ratios), we in- 
clude 8 opacities measured by iKoike et al. J2003h for the range of 0.77<x<1.0 for olivines of 
composition Mg2sFe2_2xSi04, and also x=0.21 8. In addition, we created a hybrid forsterite 
opacity using the x=l ("FolOO") opacity fro mlKoike et al.l (120031) for A < 15 fim and room- 
temperature forsterite opacity measured by IChihara et al.l (120011) for A > 15 /im. We also 
include the opacities of fo r sterit e powder measured by iFabian et al.l (1200 ll ) and annealed 
forsterite by iFabian et al.l (120001). Finally, we comp ute opacities for the Continuous Dis - 
tribution of Ellipsoids (see iBohren fc Huffman! Il983l ) and CDE2 (see IFabian et al.l I2OOII ) . 
CDE2 and CDE are distributions of ellipsoids whose size can be described as being in the 
"Rayleigh-limit" , or characteristic size (of the semimajor axes of the ellipsoid) much less than 
with ellipsoids described by every possible set of axial ratios weighted more toward near 
spherical shapes (CDE2) and equally (CDE), respectively. For both shape distributions of 
forsterite, we compute opacities usirig optic al constants for the three crystallographic axes of 
forsterite provided by ISogawa et al.l (120061 ). To account for this anisotropy of forsterite, we 
compute opacities for forsterite grains in a given shape distribution (CDE2 or CDE) three 
times, once for each of the 3 crystallographic axes' optical constants; we average these three 
resulting opacities to account for forsterite grains oriented randomly (see lBohren fc Huffman 
19831 ). 



As a starting mixture (not the same as the standard mixture determined and used 
later) of opacities for other comp onents to use i n mod eling our forsterite exemplars, we use 
the non-forsterite opacities from ISargent et al.l (120081 ): that is, CII )E2 opacities for amor- 
phous olivine and amorphous pyroxene using optical constants from iDorschner et al.l (119951 ) 
for MgFeSi04 and Mgn.7F en..'^Si0.s, respectively, "En90" ens tatite opacity (Mgn.gFeo.iSiOs) 
from IChihara et al.l (120021 ). and annealed silica opacity from IFabian et al.l (120001 ). For large 
grains, we use the aforem entioned amorphous silicate optical properties and Mie Theory 
( IBohren fc Huffmanlll983l ) to compute the opacity of 5 fim radius spheres that are 60% vac- 
uum by volume, computing effective complex dielectric functions for thi s porous amorphous 
silicate material using Bruggeman Effective Medium Theory (EMT; iBohren &: Huffman 
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19831). 



The scattering efficiencies of small amorphous and crystalline silicate grains are negligi- 
ble because they are in the Rayleigh limit for Spitzer IRS wavelengths and are not included 
in their opacities. The scattering component of opacity is not negligible for 5 /im radius 
grains. The scattering cross section is comparable to the absorption cross section in the 10 
fim region for both amorphous pyroxene and olivine dust grains. Since we model the emis- 
sion from an optically thin region (the upper disk layers), only the absorption component 
contributes to our model. Full radiative transfer treatments of the emergent intensity from 
these optically thick disks must include the effects o f scattering on t h e 5-37 /xm IRS spectra. 
We note the detailed radiativ e tran sfer models of iD'Alessio et al.l (120011 ). who cite earlier 
work by iMiyake &: Nakagawa] (Il993l ). show that scattering by large grains does not much 
change the flux from optically thin regions and makes the flux from optically thick regions 
slightly smaller. For this reason, we deem the lack of inclusion of scattering in opacities not 
of great concern. 

We give in Table 1 the mass weights in units of our conservative uncertainties, a^ (a 
measure of significance), of the olivine components for each of the seven forsterite exemplars. 
Those olivine opacities whose Sfe are highest contribute most significantly to the fit of a given 
forsterite exemplar. A given entry in a column in Table 1 is the mass weight of the component 
named in the first entry on the same row divided by its uncertainty, giving the significance 
of a particular opacity used in a given model. A prefix of "C" (or "cool" ) means the opacity 
was used at the lower temperature, and a prefix of "W" (or "warm") means the opacity was 
used at the higher temperat ure. Of the e r itries whose names include "J98" (meaning those 
opacities were presented by IJaeger et al.l (Il998l )). "fo" denotes forsterite, "ol" denotes Mg- 
rich olivine, "ho" denote s hortonolite, and "fa" denotes fayahte. The subsequent 9 opacities 
come from iKoike et al.l (120031 ). with the number after "Fo" indicating the stoichiometric 
abundance, expressed as a percentage, of Mg with respect to Fe in the ratio Mg/ (Fe+Mg) 
for olivine [Mg,Fe]2Si 04. Entries fo r FolOO, Fo90.7, and Fo77, and the opacity of forsterite 
powder presented by iFabian et al.l (120011 ) are not included at either temperature because 
no mass of any of them was ever used in any of the models of the fors terite exemplars, 
"annfor" denotes the opacity of forsterite presented by iFabian et al.l (120001 ) annealed for the 
longest duration, "fo hyb" is the hybrid op acity (explained in the text). S6 means the optical 
constants come from ISogawa et al.l ( 20061 ). and CDE and CDE2 are the shape distributions 
for which opacities are computed using these optical constants. We find that a few opacities 
recur from model to model, and some of those at both model dust temperatur es. Arrows in 



the fin al column indicate recurring opacities. Forsterite powder measured by iJaeger et al. 



( 119981 ) shows up in 4 models as cool dust, but not as warm dust. Our forsterite "hybrid" 
opacity shows up in 3 models as cool dust, but only once as warm dust. Annealed forsterite 
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occurs in 4 models as a warm dust component, but not as a cool component. Fayalite from 
Jaeger et al.l (119981 ) occurs as warm dust with frequently high a^ for the 6 exemplars, but as 
cool dust only once. The opacity that most repeatedly occurs as both warm and cool dust 
components in the models of the e xemp lars is the CDE2 opacity using optical constants for 
forsterite given by ISogawa et al.l (120061 ). Frequently associated with this is also the CDE 
opacity for forsterite from the same optical constants. 

The best fits seem to come from a c ombination of CDE2 and CDE shape distributions 
of forsterite using optical constants from ISogawa et al.l (120061 ) . suggesting that a shape dis- 
tribution intermediate in weighting between CDE2 and CDE is indicated. Further, we note 
that the CDE shape distribution is somewhat unphysical, requiring ellipsoids either infinitely 
long or infinitesimally thin. We therefore computed an opacity curve for a custom-designed 
shape distribution of ellipsoids intermediate i n weighting between that of CDE2 and of CDE. 
We note that a similar attempt was made by lZubko et al.l (119961 ) resulting in what those au- 
thors called "modified Continuous Distribution of Ellipsoids" . As the features for forsterite 
grains in the CDE shape distribution peak at wavelengths slightly longward of the peaks in 
the data and those of the CDE2 shape distribution of forsterite grains peak at wavelengths 
shortward of the peaks in the data, we expect the intermediate shape distribution to give 
opacity peaks at more optimal wavelengths than do the CDE or CDE2 shape distributions. 
Not knowing exactly how to weight one ellipsoidal shape relative to another, like CDE we 
weight all allowed shapes equally, but unlike CDE we do not allow ellipsoidal shapes with 
Lj parameter (see Bohren &: Huffman 1983 ) equal to zero in the shape distributior|^. Our 
restricted shape distribution, which we call the "truncated Continuous Distribution of El- 
lipsoids", or tCDE, weighs all points equally within the bounds of a right triangle on the 
graph bounded by the three vertices (Li, La) of (0.005, t), (0.005, 0.99), and (0.995-t, t). 
We chose not to search for the best forsterite opacity using the tCDE shape distribution of 
forsterite and an untested enstatite opacity (En90), but, instead, to search jointly for the 
best value of the bound, t, on the tCDE shape distribution and the best enstatite opacity. 
This is because enstatite opacities have features of similar widths and at many of the same 
wavelengths (10, 19, 23, 28, 33) as the opacities of forsterite. 



^The Lj parameters are inversely proportional to the semimajor axis of one of the three ellipsoidal axes, 
there is one hj parameter for each of the three axes, and the sum of the Lj parameters equals 1. Therefore, two 
Ijj parameters, Li and L2, suffice to describe completely the shape of the ellipsoid. This can be represented 
graphically with Li as the horizontal axis and L2 as the vertical axis. Any conceivable ellipsoidal shape is 
then represented on this graph as a point in the first quadrant satisfying L2 < 1-Li 
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3.3. Enstatite Opacity 



Now we focus on choosing tlie best opacit y of enstatite, also se arctiing; for the best bound 
in the tCDE shape distribution for forsterite. ISargent et aLl (120061 ) found FN Tau to require 
enstatite to fit narrow features in its 10 nm complex at 9.3, 10.6, 11.2, and 11.6 fim. In this 
study, we are interested in modeling all silicate dust features found in IRS spectra, so we 
look for guidance longward of the 10 /im complex to the spectrum of FN Tau. FN Tau has 
a very weak 16 /im feature (unlike the forsterite exemplars), a strong 28 fim feature relative 
to both the 23 and 33 /im complexes, and a double-peaked 33 fim complex. We look for 
these patterns elsewhere in our Taurus- Auriga sample to find other enstatite exemplars. DH 
Tau also has a narrow 9.3 fim feature, strong 28 /im complex, and double-peaked 33 fim 
complex. We also include Haro 6-37 and HK Tau as enstatite exemplars because both have 
narrow features at 9.3, 10.6, and also between 11 and 12 fim, though they lack prominent 
crystalline silicate features at wavelengths longward of the 10 fim complex. 

As with the olivines, numerous lab measurements of pyroxenes exist. We include opac- 
ities computed from transmission spectra of submicron pyroxene powders in potassium bro- 
mide (KBr) pellets for orthoenstatite (MgSiOs crystals in orthorhombic crystal structure), 
clinoenstatite (MgSiOs in monoclinic crystal structure), the aforementioned En90, 'En80" 
(Mgo.gFeo.aSiOg), "En70" (Mgn vFen ..Si O,), 'En60" (M gn ,Fen ,SiO.O, 'En50" (Mgo.sFeo.sSiOa), 
and "EnOO" (FeSiOs, or ferrosilite) by IChihara et al.l (120021). We also cons ider similar mea- 
surements of synthetic orthoenstatite and clinoenstatite (IKoike et al.ll2000l ): natural clinoen- 
statite from Akita, Japan, orth oenstatite from Norway, and synthetic clinoenstatite and or- 
thoenstatite ( Koike et ahlboooh: and h ypersthene (Mgo.5Feo.5Si03), bronzite (Mgo.8Feo.2Si03), 
and clinoenstatite (I Jaeger et al.lll998l ). We also compute opacities of orthoenstatite grains 
in the CDE2 shape distri b ution using optical constants for the three crystallographic axes 
provided by IJaeger et al.l (119981 ). At the same time, for the rest of the dust types in the 
models, we use the same non-forsterite and non-enstatite opacities as were used in the mod- 
els of the forsterite exemplars in the previous section, and we explore values of the bound, t, 
on the tCDE shape distribution of forsterite ranging from 0.04 to 0.1. We also include the 
opacity of the CDE2 shape distribution of forsterite grains. We use all of these opacities to 
model the spectra of the enstatite exemplars FN Tau, DH Tau, Haro 6-37, and HK Tau, in 
addition to the forsterite exemplars. 

Table 2 gives the significance of the mass weights of the opacities of crystalline pyroxene 
and crystalline olivine in the models of the enstatite and forsterite exe mplars. The opacitie s 
named "en90", "en80", "en70", "enSO", and "enOO" were presented bv lchihara et all J2002h . 
The ones named "J98CDE2" use the optical con stants for the three crystallographic axes 
of orthoenstatite presented by IJaeger et al.l (119981 ) to compute opacity for grains in a CDE2 
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shape distri b ution . "KOso" is the synthetic orthoenstatite, whose opacity is presented by 



Koike et aLl (120001 ) . "J98b" and "J98 c" are bronzite and chnoenstatite opacities, respec- 
tively, presented by iJaeger et aLl (119981 ) obtained via transmission measurements of crushed 
grain s pressed into p ellets . Entries for orthoenstatite, ch noenstatite, a nd "e n60" opacities 
from IChihara et al.l (120021). sy nthetic chnoenstatite from iKoike et al.l (120001 ). and hypers- 
thene from IJaeger et al.l (119981 ) are not included at either temperature because no mass of 
any of them was ever used in any of the models of the exemplars of both forsterite and 
enstatite modeled jo intly. "S6" in the last five rows refers to optical constants presented by 
Sogawa et al.l (120061 ) used to compute opacities for the CDE2 shape distribution (the first of 
the five) and for the tCDE shape distributions with bounding parameter (see discussion in 
Section 3.4), t, of 0.10, 0.08, 0.06, and 0.04 for the second through fifth of the five opacities. 

Among t he py roxene opacities, CDE2 orthoenstatite from optical constants provided by 
Jaeger et al.l (119981 ) shows up at significant levels as a cool dust cornpone nt but only weakly 
as a warm dust component. Ferrosilite (EnOO) by IChihara et al.l (120021 ) often contributes 
significantly as a warm dust component but not as a cool component. The one opacity that 
contributes significantly most frequently as both a warm and a cool dust co mponent is the 



chnoe nstatite (enstatite having monoclinic crystalline structure) En90 from IChihara et al. 



( I2OO2I ). Though we use the opacity of enstatite of composition Mgo.gFeo.iSiOs having slight 
iron content, the fit at long wavelengths is not ideal and suggests caution in the interpretation 



of find ing crystalline silicates with any iron content. This opacity was used bv ISargent et al. 



(I2OO6I ) in modeling the 10 /im complexes of 12 TTS and by lSargent et al.l (120081 ) in modeling 
silica exemplars. Among the various opacities of forsterite, the tCDE opacity with bound 
t=0.1 gives the best fit to both the enstatite and fors terite exempla r s. We note the improve- 
ment of fit of the spectrum of ROXs 42C reporte d bvlSargent et al.l (120081) when using tCDE 



instead of CDE forsterite (optical constants by ISogawa et al 
freedom of 4.3 and 5.1, respectively. 



20061 ). with per degree of 



3.4. Amorphous Silicate Opacities 



Now we turn to amorphous silicate opacities. ISargent et al.l (120061 ) showed that the 
transitional disks (disks whose outer parts are optically thick to mid-infrared radiation but 
whose inner regions are very optically thin) CoKu Tau/4, DM Tau, and GM Aur were fit 
well using opacities of mostly su bmicron amorphous si licates having smooth 10 and 20 fim 
features. LkCa 15 was shown by lEspaillat et al.l (120071 ) to be a pre-transitional disk, which 
are like the aforementioned transitional disks but with more dust in the inner regions. As 
with the transitional disks, LkCa 15 shows smooth 10 and 20 /xm features. FM Tau, TW 
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Cha, and UY Aur also show similarly smooth 10 and 20 fim features. We consider these 
seven objects to be amorphous silicate exemplars. 



Dorschner et al.l (Il995l ) gives complex indices of refraction for amorphous dust of py- 
roxene composition ("amorphous pyroxene") MgSiOa, Mgo.yFeo.sSiOs, and Mgo.sFeo.sSiOs 
and also for amorphous dust of olivine composition ("amorphous olivine") MgFeSi04. We 
use complex indices of refraction for samples of amorphous forsterite (Mg2Si04) and amor- 
phous enstat ite (MgSiO.':;) obtaine d by laser ablation of crystals of forsterite and enstatite, 
respectively (jScott &: Duleylll996l ). We also use complex indices of refraction for samples of 
amorphous forsterite and amorphous enstatite prepared by sp uttering of samples of MgSi 
and Mg2Si in an atmosphere of 50:50 argon-oxygen iDayl (Il979l ). Additionally, we use opti- 
cal constants for a sample of amorphous bronzite obtained by quenching a melt of a natural 



sampl e of bronzite (Mgo.gFeo.iSiOa) from Paterlestein, Germany analyzed by lDorschner et al. 



(119881 ) a nd for amorphous py roxene of cosmic composition (Can.n .'}Mgn.52Fen.45SiO.-j; sa mple 



IS") by IJaeger et al.l (119941 ). For the same reasons provided by ISargent et al.l (120061 ) . we 



assume the amorphous silicate dust grains are in the CDE2 shape distribution. 

We also compute opacities of 5 /im radius, 60% vacuum porous grains. The narrow 
features in the opacity curves of crystalline grains are very sensitive to the exact details of 
these grains' properties such as shape, composition, porosity, size, etc, as we have already 
indicated in our exploration of the opacities of submicron cr ystalline grains u s ed to fit silica, 
forsterite, and enstatite exemplars. For this reason, we, as iBouwman et al.l (120011 ) . do not 
include opacities of large crystalline grains. The opacities for large amorphous silicate grains 
are less sensitive to these details. 

Large grains (greater than 1 /xm in characteristic size) in protoplanetary disks should 
not be homogeneous. Instead, they should be like Interplanetary Dust Particles, which are 
heterogeneous aggregates whose components are ~ 0.1 /xm in size (see iHarker et al.ll2002l ) 
and cornposed chiefly of either amorphous silicate, forsterite, or enstatite (see the review by 
BradlevlEooal l How shoul d the opacities of such large heterogeneous grains appear? The 
recent study by iMin et al.l (120081 ) concluded that the opacity of a large heterogeneous ag- 
gregate is equal to the sum of opacities of its constituents, the opacities resembling those of 
homogeneous grains of characteristic size corresponding to the abundance of the particular 
constituent in the large aggregate. As an example in their study, they compute the opacity 
of a large porous aggregate of amorphous silicate, forsterite, and enstatite, the amorphous 
silicate being the most abundant component of the aggregate. The opacity of the amorphous 
component of the aggregate resembled that of homogeneous grains of amorphous silicate of 
size similar to that of the heterogeneous aggregate, while for the lowest abundance of crys- 
talline grains in the aggregate, the forsterite contribution to the aggregate opacity resembled 
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the opacity of much smaller grains of h omogeneous forsterite. Very similar conclusions were 
also reached by lBouwman et al.l (120081 ) . finding that the typical grain size of crystalline sil- 
icate grains used in their models of seven spectra were submicron, but typical amorphous 
silicate grain sizes were up to 6 fim radius (solid grains). 

Because the grains in protoplanetary disks grow primarily by sticking together, it is 
likely large grains will be heterogeneous, composed of various sub-micron components, sim- 
ilar to the IDPs. The small amorphous components are believed to come directly from the 
ISM, with sizes < 0.25 /im in radius. Since nearly all our objects show evidence for this 
amorphous component, one expects a significant amorphous component in the large, porous 
grains. Furthermore, the amorphous components were present at the initial formative stages 
of the protoplanetary disks. On the other hand, the crystalline com ponents were produced 
later, by as yet unknown processes in the protoplanetary disk s (e.g.. Ivan Boekel et al.ll2004 



Sargent et al.ll2006l : IWatson et al.ll2008l : iBouwman et al. Thus one would expect amor 

phous components to dominate the heterogeneous fluffy grains. In this case, the large grain 
spectra will resemble the sum of completely amorphous large grains with an admixture 
of sub-micron crystalline grains. The models presented here will correctly detect the large 
grains, through the broadening of the silicate features to longer wavelengths. Any crystalline 
components of large grains will be modeled as small crystalline grains. The opacity of a large 
porous heterogeneous grain with a significant abundance of crystalline components, however, 
will not be modeled well by our standard set of opacities, especially at longer wavelengths 
(A > 20 yum) in the IRS spectra (see Appendix B). Because of optical depth effects, this will 
somewhat underestimate the crystalline mass fraction contributed by large grains. 

To model our amorphous silicate exemplars, we use annealed silica, enstatite, and 
forsterite in the tCDE shape distribution (with bound t=0.1) as our crystalline silicate 
opacities. We use the aforementioned amorphous silicate optical properties for grains in the 
CDE2 shape distribution to compute opacities of submicron grains of amorphous pyroxene 
and amorphous olivine and also of 5 fim radius, 60% vacuum grains of amorphous pyroxene 
and amorphous olivine. In testing the amorphous pyroxene and amorphous olivine opacities, 
we chose the material whose opacities of both submicron and 5 fim grains showed up most 
frequently as a significant contributor as both warm and cool dust. In Table 3, we list the 
significance of the mass weights of the amorphous pyroxenes and amorphous olivines used 
in the models of our seven arn o rphou s silicate exemplars. "D95" refers to optical constants 
presen ted by iDorschner et al.l (119951 ) . "D79" refers to o ptical constants prese nted by iDay 
Jl979h . "SD96" refers to optical constar its presented bv IScott fc DulevI Jl996h . and "J94" 
refers to optical constants p resented by IJaeger et al.l (119941 ). The amorphou s olivine "0 1" 
from IDorschner et al.l (119951 ) is MgFeSi04, and the amorphous olivine from iDayl (119791 ) is 
amorphous forsterite. The amorphous pyroxenes with "Py5" and "PylO" in their names 
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refer, respectively, to Mgn 5Fen.5SiO.ci and MgS iO.s presented by lDorschner et al.l (Il995l ). The 
amorphous pyroxene from I Jaeger et aL is amor phous pyroxene of cosm ic composition, 

and the amorphous pyroxenes from iDayl ( 19791 ) and IScott &: Duleyl ( 19961 ) are amorphous 
enstatite. "Sm" denotes small Rayleigh-limit (here, submicron) size grains in the CDE2 
shape distribution. "Lg" denotes large 5 /im radius 60% vacuum porous s pheres. Entries 
for op acities computed from optical constants for amorphous forst erite fromlScott fc Dulev 



(Il996l ). amorphous pyroxene of composition Mgn.yF en .sSiO.'H from iDorschner et al.l (119951 ). 



and amorphous bronzite from IDorschner et al.l (119881 ) are not included at either temperature 
because no mass of any of them was ever used in any of the models of the amorphous silicate 
exemplars. 

Although it is ne v er use d as a cool dust component, the amorphous olivine MgFeSi04 
from IDorschner et al.l (|l995l ) is the best amorphous olivine because it is a very significant 
warm dust component in the models of all seven amorphous exemplars. LkCa 15 provides 
one of the most stringent tests of the amorphous silicate feat ures, as its silicate feat ures are 
very prominent and have high S/N ratios. Amorphous olivine (IDorschner et al.lll995l ) almost 
exclusively fits this spectrum's 10 fim feature. Amorphous pyroxene can dominate as a cool 
dust compon ent over amorphou s olivine; for instance, the amorphous pyroxene of cosmic 
composition (jjaeger et al.lll994j ) exclusively fits the 20 /xm feature of LkCa 15. Although 
the rest of the amorphous pyroxenes are fairly equal conten ders for best a r norph ous pyroxene 
opacity, the amorphous pyroxene of cosmic composition by IJaeger et al.l (119941 ) shows up as 
the most significant amorphous pyroxene opacity in the most exemplars, so we choose this 
as our amorphous pyroxene opacity. The minor modification to the amorphous pyroxene 
opacity usually improved (but sometimes made wors e) the by ~ 0.2 over the amorphous 
pyroxene opacity used by ISargent et al. J2006l . l2008h . 



The amorphous pyroxene of cosmic composition used in our models has a magnesium-to- 
iron ratio of 52:45, while the amorphous olivine used has a Mg-to-Fe ratio of 50:50. These are 
not far from the Mg -to-Fe ratio for the " c osmic" abundances of these tw o elements (Mg:Fe 
~ 4:5, measured bv lHolweger et al.lll990l : ISnow fc Wittlll995l : ICoxIboOOl ). We also note the 
imaginary parts of the complex dielectric functions in the visible and near-iri frared regions 
(0.2-8 fim wavelength s) of the best amo rphous olivine (IDorschner et al.lll995l ) and the best 
amorphous pyroxene (I Jaeger et al.lll994j ) are very similar to that of "astronomical silicate" 
(IDraine fc Leelll984l ). whose near-ultraviolet, visible, and near- infrared imaginary part of 
the dielectric function was specifically constructed so that astronomical silicate grains would 
heat to the correct temperatures in the presence of stars. This was important in attempting 



to obt ain a self-consistent model of the Spitzer IRS spectrum of IP Tau by ISargent et al. 
(I2OO6I ). We use the best opacities that fit our various forsterite, enstatite, and amorphous 



silicate exemplars to fit our sample of 65 spectra of TTS from the Taurus-Auriga star-forming 



- 18 - 



region. 



The first real test of tliis ensemble of dust opacities to be used as a standard mixture for 
dust emission models was to fit the mi d-infrared spectrurn of the Trape zium nebulosity. The 
8-13 fim spectrum of the Trapezium (IForrest et al.lll975l ) was used by lDraine fc Led (119841 ) 
to derive the emissivity of "astronomical silicate" grains, which represent th e dust in the 



Inters tellar Medium. The 8-13 fim spectrum of the Trapezium presented by iForrest et al. 



(119751 ) was combined with the 16-38 /xm spectrum of the Trapezium obtained by subtract- 
ing 3% of the spectrur n of the Kleinman- Low nebula from the spectrum obtained pointed 
at the Trapezium (see IForrest et al.lll976l ). This accounted for contamination by emission 
from the KL nebula of the wide beam of the detector while pointed toward the Trapezium 
(IForrest et al.lll976l ). The best-fit model for the Trapezium spectrum is shown in Figure 1. 



As can be seen, the dominant components used to fit the spectrum are submicron 
amorphous pyroxene and submicron amorphous olivine. This is consistent with more recent 
analyses of dust com position of the ISM along the line-of-sight to the Galactic Center by 
Kemper et al.l (120041 ). who found their dust models required negligible large grains and very 
small amounts of crystalline silicates (less than 2.2% by mass). The parameters for the model 
shown in Figure 1 are given in Table 4. 

We describe in the Appendix our tests of Bruggeman effective medium theory used to 
compute effective complex dielectric functions and Mie Theory used to compute opacities of 
large grains. 



3.5. Degeneracy of Model Components 

We desire to measure the degeneracy between model components. This is the extent to 
which one component could be replaced by other components and achieve a similarly good 
model fit. After computing the best-fit models (to be described in the next chapter) for the 
spectra of FN Tau (an enstatite exemplar), IS Tau (a forsterite exemplar), ZZ Tau (a silica 
exemplar), and DM Tau (an amorphous silicate exemplar), was minimized with respect 
to the blackbody solid angles and dust mass weights over 7.7-37 fim wavelengths at the 
two dust temperatures (found by the best-fit models of these exemplars), not eliminating 
components with negative solid angle or mass weight. In the process of minimizing x^, a 
16x16 element covariance matrix was computed, each of the rows and columns of which 
belong to one of the 7 dust species or the blackbody at one of two temperatures. Each of the 
elements of this matrix is a sum over all concerned wavelengths (7.7-37 fim) of a product. 
The product is opacity times the Planck function (or just the Planck function in the case of 
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a blackbody component) specific for a given row times tlie opacity times tlie Planck function 
(or, again, just the Planck function in tlie case of a blackbody component) specific for a 
given column, all divided by the square of the flux density uncertai nty. The inverse of this 



matrix gives the covariance matrix (see discussion in Chapter 14 of iPress et al.lll992l ). The 
diagonal elements of this matrix are the variances of each of the dust components in the 
model, and the off-diagonal elements are the covariances of the dust component of a given 
row with the dust component of a given column. 

The correlation coefficient, r, is computed for each off-diagonal element by dividing the 
covariance of that element by the square root of the product of the variances of the two dust 
components in question. One variance is the diagonal element of the same column, the other 
is that of the same row. Two highly degenerate components will be highly anticorrelated, 
with a correlation coefficient very near -1. A correlation coefficient near zero means the two 
components are not correlated. In addition, for each off-diagonal element, the probability, P, 
of a coefficient of equal or greater magnitude being found for a non-correlated data set (which 
is t he probability of the correlation coefficient having been drawn from a random distribution; 



see lTaylorlll982l ) is computed. P near 0% indicate significant correlation coefficients, and P 
near 100% indicate insignificant correlation coefficients. 

Here we discuss the most significantly degenerate component pairs, which are shown in 
Table 5. The most negative r's are for the pair of cool large amorphous olivine and cool 
large amorphous pyroxene. For this pair, r is between -0.87 and -0.89 (and P=0.0%) for the 
four representative exemplars FN Tau (enstatite exemplar), IS Tau (forsterite exemplar), 
ZZ Tau (silica exemplar), and DM Tau (amorphous silicate exemplar), indicating the most 
degeneracy between components. Warm large amorphous olivine and warm large amorphous 
pyroxene are also fairly highly anticorrelated (degenerate); for this pair, r is usually around - 
0.72 over the sample of 65 spectra. Similarly, cool small amorphous pyroxene and cool small 
amorphous olivine are fairly significantly degenerate with each other, with -0.82>r>-0.89 
over the four representative exemplars (P=0.0%). The warm small amorphous pyroxene and 
amorphous olivine are also highly degenerate, with r usually near -0.73 over the sample of 65 
spectra. Cool enstatite is typically fairly significantly degenerate with cool forsterite, with 
r usually near r=-0.46 over the sample of 65 spectra, as enstatite and forsterite share very 
similar features in their opacities for wavelengths longward of the 10 /xm complex. Silica 
also shares a 20 fim feature close to 19 fim complex of enstatite, explaining the spread of 
-0.29>r>-0.50 (P<0.1%) in this pair of components for the four representative exemplars. 
At warm temperatures, enstatite is only significantly degenerate with silica, with r for the 
65 spectra in the sample usually being about -0.50, as warm silica and warm enstatite 
have strong, narrow features at 9.3 and 9.1 /im, respectively, while warm forsterite is not 
degenerate with either, not sharing a strong feature with silica or enstatite in the 10 /xm 
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region. 

For the highly degenerate pairs, the dust emission models accurately determine the sum 
of the masses of the two components, but not so accurately the individual masses. This 
indicates the model finding that the inner disk (warm component) tends to be dominated by 
amorphous olivine while the outer disk (cool component) shows more amorphous pyroxene is 
probably not physically real. This, at least partially, explains the negative masses often found 
to give the very best fit. By zeroing the negative mass, the complementary component's mass 
is increased to give nearly as good a fit. Figure 2 demonstrates this effect by showing the 
degeneracy between cool small amorphous pyroxene and cool small amorphous olivine for 
the model of the Trapezium shown in Figure 1. For the best-fit model temperature pair, 
the model is recomputed but without setting any of the components in the model to zero, 
allowing mass weights to be negative if that results in the lowest per degree of freedom 
(reduced x^)- All resulting model parameters are held at their new best-fit values except for 
the mass weights of cool small amorphous pyroxene and cool small amorphous olivine, which 
are varied over the ranges indicated on the plot. The levels of reduced resulting from 
exploring this range of mass weights of these two dust components are shown as contours 
on the plot, with the levels of the first three contours indicated. This contour plot suggests 
that replacing some amount of one of the cool small amorphous silicates with the same mass 
of the other gives a very similarly good fit, changing reduced only very slightly. 

Histograms of the correlation coefficient between four representative pairs of dust com- 
ponents for all 65 spectra in the sample arc given in Figure 3. The histograms all show single, 
well-defined peaks. This indicates the correlation coefficients are measuring real degeneracy 
between model components. 

4. Results 

Dust emission models are fit to the 65 Taurus-Auriga spectra, which we show in Figures 
4-14. Table 6 gives the parameters of these models (temperatures, blackbody solid angles, 
mass weights, reduced x^). x^ is minimized over 7.7-37 /xm, so both dust temperatures used 
in the models are well constrained by relative uncertainties of 10%. The median high and 
low model dust temperatures are 545 K and 127 K, respectively. In the DR Tau model, which 
has dust at precisely these temperatures, most of the 10.0 //m wavelength flux above that 
from the blackbodies originates from 545 K dust, while most of the 20.0 //m flux above that 
from the blackbodies originates from 127 K dust. Assuming a distance of 140pc to DR Tau, 
assuming all of the optically thin dust fiux at 10.0 /xm comes from 545 K dust in a disk of 
total optical depth at 10.0 \im. of 0.1, and assuming all of the optically thin dust flux at 
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20.0 /iin comes from 127 K dust in a disk of total optical depth at 20.0 /im of 0.1, the radii 
of these two disks in the case of DR Tau are 0.75 AU and 11.5 AU for the warm and 
cool optically thin dust, respectively. We list at the end of Table 6 the best opacities that 
fit our various forsterite, enstatite, and amorphous silicate exemplars to fit our sample of 65 
spectra of TTS from the Taurus- Auriga star-forming region. 



4.1. Extreme Inner Disk Grain Growth 



The spectra of DM Tau, DO Tau, UZ Tau/e, XZ Tau, and ZZ Tau IRS in our sample 
show evidence of high abundances of large grains. Figure 4 shows our models of these 
five spectra, breaking into components the model of UZ Tau/e as an example of the models. 
Warm large amorphous olivine dominates the 10 /im complexes of these five objects' spectral 
models. The 20 /xm features also are partly fit by this same component, but usually more of 
the fiux required to fit the 20 /im features comes from cool submicron amorphous pyroxene. 
The overall shape of the spectrum of DM Tau differs f rom that of the other spectra in Figure 
4 because it is a transitional disk f Calvet et al.|[2005l ). lacking significant dust closer to the 
central star than ~ 3 AU. We note that these five spectra support our use of opacities of 
large grains of amorphous silicates and do not require those of large crystalline grains. That 
our models fit our spectra well is consistent with amorphous silicates being more abundant 
both as small grains and as constituents in larger grains (see Section 3.4). We note that 
significant abundanc es of large grains hav e been fou nd for SR20, a TTS in the Ophiuchus 
star-forming re gion ( AfcClure et al.l [ 2008h CS Cha jEspaillat et al.lboOTh . and CVSO 224 
Espaillat et al.l (120081 ) . iFurlan et al.l (120071 ) found a large abundance of large grains around 
HD 98800 B, a somewhat older (~ 10 Myr old) YSO in the TW Hydrae as sociation. Large 



abund ance s of large grains around othe r somewhat older stars were found by lBouwman et al. 
(l2008h and lKessler-Silacci et all J2006h . 



4.2. Prominent Forsterite Spectral Features 

In Figure 5 we show the spectra of six spectra with prominent forsterite features, five of 
which are our forsterite exemplars, with the model of IS Tau broken down into components. 
We note in support of the forsterite opacities used in our models that the fits at all wave- 
lengths to the spectra of these stars, DK Tau, F04147+2822, GN Tau, IS Tau, V836 Tau, 
and V955 Tau are quite good. All major features and complexes of forsterite at 10.0, 11.2, 
16, 19, 23, and 33 /im are fairly well fit. 
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The spectrum of F04147+2822 (Figure 5) is especially interesting. Its 10 /im complex 
lacks any of the distinctive narrow emission features characteristic of the forsterite, enstatite, 
or silica and is fit well by a combination of emission from submicron amorphous silicate grains 
and large amorphous pyroxene. At longer wavelengths, however, emission features at 19, 23, 
and 33 ixm characteristic of forsterite dominate and, correspondingly, are fit well by our 
forsterite profile for grains in a tCDE shape distribution. This suggests a greater abundance 
of forsterite in the outer disk than the inner disk. One possible explanation is that this 
spectrum is a sum of emission from two disks, one with very little crystallinity throughout 
and blue in continuum color and another disk with a lot of crystallinity throughout the disk 
and red in continuum color, but this system is not known to be multiple. I f this spectrum is 



from only one protoplanetary disk, this would contrast with the finding by Ivan Boekel et al. 



(120041 ) of greater crystallinity within 2 AU than outside of 2 AU in disks around 3 Herbig 
Ae/Be stars. However, the two are not inconsistent, as T Tauri stars are less luminous than 
Herbig Ae/Be stars. Two AU in HAeBe stars should correspond to a much lesser radius 
in disks around T Tauri stars. It suggests that, on average, the regions to which we refer 
as "inner" and "outer" disk regions in the population of 65 T Tauri stars in Taurus- Auriga 
are both outside of the regions in T Tauri star disks analogous to the region inside 2 AU of 
Herbig Ae/Be disks. 

In some spectra, however, the long-wavelength side of the 10 /im feature is not fit well. 
DD Tau, DE Tau, and V710 Tau all share this problem, as is seen in Figure 6, with the model 
of FX Tau broken into its components. There are many ways by which forsterite opacity 
features can be centered at longer wavelengths than those in our tCDE opacity curve. Grain 
shape distributions more heavily weighted toward extreme shapes (extremely fiat, extremely 
elongated), larger grains (or large aggregates of small grains, with a large fraction of these 
small grains being forsterite), olivine grains with nonzero Fe/Mg ratios, and forsterite grains 
with greater porosity would all result in forsterite emission features at longer wavelengths 
than resulted from the forsterite tCDE curve (solid spheres of forsterite give rise to opacity 
features at shorter wavelengths than those of forsterite in the CDE shape distribution, and 
the opacity curve of 60% po rous spheres of forste rite is almost identical to that of forsterite in 



the CDE shape distribution ISargent et al.ll2006l ). Also, both CoKu Tau/3 and FX Tau (Fig 



4) have 23 and 33 /im emission complexes with emission that extends to longer wavelengths 
than provided by the submicron forsterite grains in the tCDE shape distribution in our 
models. That both spectra lack significant 28 /im complexes indicates the dominant form of 
crystals giving rise to 20-37 /im emission for these two systems is forsterite and not enstatite. 
It should be noted, though, that the 10 /im complexes in the spectra of both CoKu Tau/3 
and FX Tau are fit fairly well. We therefore conclude there is variation in the exact details of 
the olivine grain populations in the protoplanetary disks in our sample, in terms of average 
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grain shape, grain size, grain composition, and grain porosity. For DD Tau, DE Tau, and 
V710 Tau, this variation occurs for the warmer forsterite grains, while for CoKu Tau/3 and 
FX Tau, this variation occurs for the cooler forsterite grains. We note that most of these five 
also seem to have emission features at 14 /im wavelength (see below. Section 4.8). Whatever 
gives rise to these 14 /xm features is not necessarily related to the poor fits to the forsterite 
features because the poor fitting happens over very limited wavelength ranges of 11-11.5 /im 
and 23-25 /xm. 



4.3. Variation in Silica 



In addition to the olivine component of our models, the silica component may also vary 
in details of grain properties among all disks in our sample. While the silica features at 9, 
12.6, and 20 /im in the spectrum of ZZ Tau (and the other silica exemplars) are fairly fit well 
with the annealed silica opacity used, others like the 12.6 /im features in the spectra of DN 
Tau and FZ Tau are not quite so well fit (see Figure 7). FZ Tau has a single-peaked feature 
in this wavelength range, but it peaks at 12.45 fim, shortward of the wavelength at which 
annealed silica peaks. This could suggest the polymorph of silica giving rise to this feature is 



g-qu artz at somewhat elevated temperatures of 500-625 K (see discussion by ISargent et al. 



20081 ). D N Tau appears to ha ve a double-peaked feature, which would suggest a-quartz at 
~ 300 K (ISargent et al.ll2008l ). but the S/N at these wavelengths is low for this spectrum 
so its reality is more doubtful (not e the discussi o n on the effect of unresolved emission 
lines in this section). According to ISargent et al.l (120081 ) . the presence of silica argues for 
high-temperature processing, as it is not present in the interstellar medium and therefore 
should not be present as the starting dust mixture as the protoplanetary disk formed from 
its protostellar envelope. If a-quartz can be confirmed, the silica, once formed, must cool 
slowly enough to allow phase transition from a higher-temperature polymorph, like tridymite 
or cristobalite, to /5-quartz, then a-quartz. 



4.4. Enstatite Exemplars 

We show the model fits to our enstatite exemplars in Figure 8, with the model of FN Tau 
broken into components. Other than a mismatch between model and spectrum for FN Tau at 
11.6 fim, the fits to the 10 fim complexes are generally good. However, we note the poorness 
of fit of our models to the longest wavelength complexes (28 and 33 fim) of the spectra of 
FN Tau and DH Tau, our two best enstatite exemplars. The height-to-continuum ratios in 
our models of these two exemplars for the 28 and 33 fim complexes are unacceptably low 
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compared to those in the spectra. We note the progressively increasing height-to-continuum 
ratio of the 19, 23, 28, and 33 fim complexes in the spectra of both DH Tau and FN Tau. 
The rest of our sample are sufficiently fit using crystalline dust at two temperatures. This 
suggests that a population of submicron enstatite or forsterite dust at a third (very low) 
temperature would be required to give acceptable fits to these spectra. We note this lack of 
a third model temperature is not a problem for spectra requiring abundant forsterite (Figures 
5 and 6). Note in Figure 6 the problem is not insufficient peak-to-continuum ratio for the 
model forsterite features; rather, the problem is insufficient width of the features, which is a 
problem with the forsterite opacity. 

The LH part of the spectrum of Haro 6-37 (19.3-37 /^.m wavelength) suffers from an 
artifact, in which each of the 10 spectral orders are tilted such that the flux density of their 
short- wavelength end is higher than it should be and the flux density of their long- wavelength 
end is lower than it should be. These tilted spectral orders, which could be interpreted to 
resemble narrow crystalline silicate features longward of 19.3 fim wavelength, result from an 
artifact of the reduction of the LH spectrum of this object. This artifact could be due either 
to unsubtractcd sky emission or to mispointing of the telescope with respect to Haro 6-37 
along or perpendicular to the length of the LH slit. 



4.5. Amorphous Silicate Exemplars and Transitional Disks 

Spectra requiring abundant submicron amorphous silicate grains and model flts to these 
spectra are shown in Figures 9 and 10. Figure 9 shows the spectra and corresponding models 
of five of our amorphous silicate exemplars (see Section 3.4) and the model breakdown into 
components for FM Tau. In general, from 7.7 to slightly longward of 30 fira wavelength, 
the fits of the models to the data are excellent. For GM Aur and especially LkCa 15, the 
each model flts the continuum at wavelengths greater than 30 /im with emission from a 
weak 33 /im complex of forsterite. This is an artifact, as these two spectra require a third 
model temperature to flt the data using emission from a blackbody or amorphous silicates 
at a very low temperature. Figure 10 shows spectra whose models also require relatively 
high abundances of submicron amorphous silicates, though they require more crystalline 
silicates and large grains than those shown in Figure 9; the model of HQ Tau is broken into 
its components as an example. Again, the model fits to the data are very good, with the 
exception of a slight insufficiency in the peak-to-continuum ratio of the forsterite 33 /xm 
complex in the model of DP Tau. 
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4.6. Small-Equivalent- Width 10 /iin Complexes 

The spectra of CY Tau, DF Tau, DG Tau, DL Tau, DQ Tau, IT Tau, and V807 Tau show 
10 /iin complexes of small equivalent width (see Figure 11; the model of DF Tau is broken into 
components). CY Tau, DF Tau, DL Tau, and IT Tau also have essentially fiat spect ra over all 



IRS wavelengths, being consistent with the models presented by lFurlan et al.l (120061 ) and their 
interpretation of settling of dust grains. Settling of dust in a disk results in the disk being 
flatter and less flared, which gives rise to a bluer spectral contin uum color in the mid-infrared 



and also less equivalent width in the distinctive silicate features. [Andrews fc Williams! (120051 ) 
noted CY Tau for its fiat spectral colors over infrared wavelengths. 

These spectra are well fit mostly by emission from submicron dust grains and the two 
blackbody components in the models. CY Tau requires a small amount of silica to fit in the 
10 fim complex and a moderate amount of cool forsterite to fit longer wavelength complexes, 
especially the 33 fim feature. The emission required by large grains in the model for DF Tau 
to fit the longer wavelengths of its spectrum is the largest of the five low 10 fim equivalent 
width spectra, in addition to small amounts of emission from submicron dust grains. Little 
emission from forsterite is required to give rise to the 10 /im feature of DG Tau, while a fair 
amount of submicron amorphous pyroxene is required to fit its mild 20 fim feature. The 
10 /im feature of DL Tau requires modest amounts of emission from forsterite and silica 
and lesser amounts from other kinds of dust. The LH spectrum (19.3-37 /im) of DL Tau 
suffers from an order-tilt artifact as does Haro 6-37 (see above), limiting the precision with 
which we can determine the crystallinity and large grain content of the cooler outer disk, 
but it appears to be optimally fit using a mixture of both amorphous and crystalline grains. 
DQ Tau has subtle features at 9.3, 9.8, 10.6, 11.1, and 11.6 /im that require a moderate 
amount of enstatite. IT Tau requires small amounts of many of the submicron amorphous 
and crystalline grains to fit its miniscule 10 /im complex, but its broad 20 /im feature is fit 
well by submicron amorphous pyroxene. The same applies to the fit of the 20 /im feature 
of V807 Tau, though it requires a modest amount of large amorphous pyroxene to fit its 10 
/im feature. The abundance of crystalline grains with respe ct to amorphous gra ins in the 



models of these five spectra is consistent with the finding by IWatson et al.l (120081 ) using all 
(and more) of the spectra analyzed in this study that increased crystallinity accompanies 
more advanced settling of dust in disks as measured by continuum indices ne-ia and ni3_3i. 
It is also noted that these seven spectra with small equivalent width 10 /im features do not 
indicate large abundances of large grains (see Table 6). We explore these issues more in 
Section 5. 
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4.7. Mixed Dust Compositions 

In Figures 12-14, we show the spectra and corresponding models to the rest of our 
sample, with eight spectra in each figure arranged top-to-bottom alphabetically by their 
names. For all of these spectra, the abundances of the various dust species required by the 
models are mixed. Typically, no one dust type has an abundance required by its model much 
more than the abundances of any of the other dust types in that model. The spectrum of 
AA Tau (and others) have what appe ar to be a broa.d emi ssion feature at 14 /im, but this 



probably originates from gas emission (ICarr &: Najitall2008l ): also, SH and LH spectra of AA 



Tau and a few other objects show unresolved emission lines from H2O and OH (see next 
subsection). The models of FP Tau, GI Tau, and VY Tau are broken into components in 
Figures 12, 13, and 14, respectively. Note that FS Tau and FV Tau in Figure 13 and GK 
Tau, HN Tau, and Haro 6-28 in Figure 14 suffer from the same LH order-tilt artifact as did 
Haro 6-37 (see Section 4.4). Also note the high point-to-point noise in the 20-25 fim parts 
of the LL spectra of FT Tau, GG Tau B, and GH Tau in Figure 13 and IQ Tau in Figure 
14. This "fringing" artifact is due to the delamination of the LL order-sorting filter and 
shows up when the telescope has been mispointed with respect to a source in the direction 
perpendicular to the slit for a LL observation, with increasing severity of the artifact for 
greater mispointing of the telescope. 



4.8. Unresolved Emission lines 



There are some very narrow gaseous emission lines in our spectra. The spectra of 
DG Tau (Figure 11), DM Tau (Figure 4), and FS Tau (Figure 13) all show 12.8 /im fea- 
tures arising f r om [N ell]. This can affect identificatio n of the silica compo nent, as noted by 
Sargent et al. ( 20081 ) with respect to the modeling by Sargent et al. (2006) of TW Hya. The 
presence of t his feature in additio n to a narrow feature at ~ 12.4 /im in the SL spectrum of 
TW Hya led lSargent et al.l (120061 ) to conclude the presence of a-quartz, which has a double- 
peaked feature matching these two wavelengths. Investigation of high-resolution spectra of 
TW Hya showed that these two features belonged not to a-quartz but to the Humphreys-a 
Hydrogen line (HI n=7-6) at 12.37 /im and [Nell] at 12.81 yum. 

C2H2 appears in the spectra of CoKu Tau/3 (Figure 6), DF Tau (Figure 11), and DL 
Tau (Figure 11). Thi s emission line is located at the more innocuous wavelength of 13.7 fim 
(ICarr fc Najitall2008l ). and does not appear to have significantly affected the dust models of 
the spect ra of these three obj ects. However, emission from HCN centered just longward of 
14.0 yum (jCarr fc Najitall2008l ) is more problematic. Examples of spectra with this feature are 
AA Tau (Figure 12), BP Tau (Figure 12), and IT Tau (Figure 11). First, the feature occurs 
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often in spectra taken with the SL and LL modules, which we sphce together at precisely 
14.00 fim. Often there is a flux level mismatch between SL and LL at this wavelength, 
making establishing the reality of this feature difficult. Second, the width of this feature 
varies. We have no dust component with an emission feature centered around 14 /xm, so 
when the feature is wider, more data deviate from our model, worsening t he model fit to the 



data. Thirdly, observations of evolved stars in the Large Magellanic Cloud (ISloan et al.ll2006 
20081 ) sometimes show a feature at ~ 14 /im probably belonging to dust related to silicates. 



Lastly, there are emission features near the se wavelengths in spectra of Calcium-Aluminum 
Inclusi ons in meteorites (jPosch et al.l 120071 ) and in laboratory spectra of the melilite solid 
series (jChihara et al.l 120071 ). making identification of all of the 14 /im features with HCN 
emission more ambiguous. 

Finally, emission lines of H2O have been d iscovered in high reso lution Spitzer IRS spectra 
of classical TTSs in Taurus, such as AA Tau (jCarr fc Naiitall2008l ). These lines are found in 
its SH and LH spectra, especially past 25 um. Many of these lines were found iri the S pitzer 
IRS spectrum of the Class YSO protostar NGC 1333 IRAS 4B JWatson et al.ll2007h . The 
lines are generally spread out over wavelength, but there are clumps of lines around 24.5-25.5, 
27.5-28.5, and 30-31 /im wavelength. These water lines are unresolved, and our SH and LH 
spectra (R ~ 600) have been rebinned to SL/LL resolution (R ~ 90), so the effect of the 
lines generally being distributed over all wavelengths in SH and LH is to raise slightly the 
level of the continuum above the continuum underneath the water lines in the high resolution 
spectra. The clumps at 25, 28, and 30.5 fim in the original high-resolution spectra, when 
rebinned, will show up as very slight ~ 0.5 fim wide bumps in the rebinned spectra. We 
typically see no prominent crystalline silicate features centered at either 25 or 30.5 /xm in 
our spectra (though DH Tau has a small f eature centered arou nd 30.5 /im wavelength that 
has been attributed to enstatite before; see lMolster et al.ll2002l ). but we do see many 28 fj,m 
features in our data. Both the opacities of forsterite and enstatite have 28 fim complexes, 
the 28 fim complex of enstatite being more prominent with respect to its 23 and 33 fim 
complexes, so the net effect of not accounting for the water emission might be obtaining a 
enstatite abundance slightly too large. 



4.9. Imperfect Extinction Correction 

Though we attempted to correct for extinction as explained in Section 2, two spectra 
proved to be difficult in this regard, those of FV Tau (Figure 13) and CoKu Tau/3 (Figure 
6). For FV Tau, Av=5.33, and for CoKu Tau/3, Av=5. Both were just under the upper 
limit of Ay =6 above which we would not correct for extinction and therefore not model, and 
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both spectra show absorption at 15.2 fim from CO2 ice; additionally, CoKu Tau/3 has ice 
absorption feat ures at 6.0 and 6.8 A t m wavelengths as do heavily embedded C lass I YSOs in 
Taurus- Auriga JWatson et al.lbood : IZasowski et"aDl2007l : ISoogert et aDl2008h . Though the 
model fit to the 10 /im complex of CoKu Tau/3 is pretty good, the same is not true for FV 
Tau. We suspect this to result from imperfect extinction correction. At longer wavelengths, 
however, the fit to FV Tau (excluding the region around the 15.2 /im CO2 ice absorption 
feature) is quite good, while the fit to the crystalline silicate complexes at 23 and 33 /xm is 
only adequate. As any extinction correction at longer wavelengths is less than that for the 
10 /im silicate complex, we attribute this not to a problem with extinction correction, but 
rather to the opacity used to fit forsterite, which is discussed in the next subsection. 



5. Discussion 

5.1. Inner versus Outer Disk Crystallinity and Grain Growth 

From the models for all 65 spectra, we computed the percentage of mass in a given 
dust grain species at one temperature out of all mass in dust at that temperature. With 
these mass percentages, we compute histograms displaying both warm and cool dust mass 
percentages. The histograms for total warm and cool large grain mass fraction, total warm 
and cool crystalline grain mass fraction, warm and cool forsterite mass fraction, warm and 
cool enstatite mass fraction, and warm and cool silica mass fraction are given in Fig. 13. 

The test described and shown in Appendix B of fitting a 5 /im radius heterogeneous 
grain profile with our standard model suggests that the abundance of crystalline grains could 
be underestimated by our standard models by up to 50%. In addition, it also suggests the 
large grain abundance could be slightly underestimated by 23%. The net effect of this bias 
is to underestimate the true crystalline and large grain abundances. 

Figure 15a shows that the inner disk regions typically have much larger mass fraction 
in large grains than the outer disk regions. The mean and median in the histogram for 
warm large grains are around 50%; for cool large grains, they are between and 10 percent, 
though with a considerable "tail" up to 90%. Noting t hat our 5 /im radius porous grains are 



very similar in shape of opacity to 2 /im solid grains (jSargent et al.ll2006l ). our average and 



median warm large grain mass fract i ons o f 44% and 45%, respectively, compare favorably to 



the finding by ISicilia-Aguilar et al.l (120071 ) of the average grain size for 1-2 Myr old systems 
being between 1.5 and 2 /im radius (as probed by the 10 fim feature). Our average and 
median cool large grain mass fraction were 17% and 0%, respectively. 



Figure 15b shows histograms of warm and cool crystalline dust mass fraction. Both 
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histograms peak at low percentages of about 5% for warm crystalline dust and about 12% 
for cool crystalline dust (each of these percentages being the mode of the distribution), 
and there is considerable overlap. The cool crystalline dust histogram extends slightly to 
higher mass fraction bins, but we note the greater large grain mass fractions for the warmer 
inner disk regions. The average and median mass fractions for warm crystalline dust were 
17% and 11%, respectively, and the average and median mass fractions for cool crystalline 
dust were 23% and 15%, respectively. Greater large grain mass fractions can lead to slight 
underestimation of the crystalline abundance, so the amount of crystallinity in the inner 
disk regions could be higher and closer to that of outer disk regions. Our findings of the 



mean, median, and mode of the warm crystal 
5% are consis tent with those by 



ionda et al 



i ne du st mass fractions of 17%, 11%, and ~ 
(I20061). but th ey are slightly less consistent 



with those by ISicilia-Aguilar et al.l (l2007l ). iHonda et al.l (l2006l ) find typical warm crystalline 
dust mass fractions (as probe d by the 10 /xm feature) ty pically to be between 5 and 20%, 
regardless of system age, while ISicilia-Aguilar et al.l (120071 ) find typical warm crystalline dust 
mass fractions to be between 5 and 7%. Though crystalline mass fractions may be greater 
in the inner disk regions than outer disk regions (see Table 6), the total mass of the cool 
dust in our models is between ~ 10 and 1000 times greater than that of the warm dust in 
the models. The total mass in crystalline silicates is often greater for the cool crystalline 
dust than the warm crystalline dust, as Table 6 shows for V955 Tau. The crystalline mass 
fractions for warm and cool crystals are 24.7% and 20.6%, respectively, being very similar. 
However, there is ~ 250 times more mass is cool dust than warm, so there is much more 
mass in 115 K crystalline dust than in 488 K crystalline dust according to the model. 

Figures 15c, 15d, and 15e show the enstatite, forsterite, and silica mass fraction his- 
tograms. In Figure 15c, the warm enstatite mass fractions are very slightly displaced to 
larger values than the cool enstatite mass fractions. The previously discussed potential for 
slight underestimation of crystalline grain mass fraction for the inner disks due to greater 
grain growth reinforces our conclusion that the abundance of enstatite is greater in the in- 
ner disk regions. Figure 15d shows that both the warm and cool forsterite mass fraction 
histograms show a peak in the lowest bin; however, a large fraction of the cool forsterite 
distribution extends to greater mass fractions. First, we note the larger uncertainties in cool 
forsterite mass fractions than warm forsterite mass fractions (Table 6). Second, we once 
again note the potential to underestimate inner disk crystallinity. Together, these two con- 
siderations suggest the warm forsterite abundance is similar to that of cool forsterite for our 
sample. Figure 15e shows the same is true for silica as for forsterite; for the same reasons as 
forsterite, the warm and cool silica abundances in our sample are similar. The higher relative 
abundance of enstatite over forsterit e in inner disk regions than for outer disk regions was 
also found by iBouwman et al.l (120081 ) . 
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5.2. Correlations 



In order to search for correlations between pairs of model parameters, stellar properties 
and disk properties, we have computed the linear correlation coefficient, r ( BevingtonI 1969 ) 



and its corresponding probability, P, of finding a linear correlation coefficient of magnitude 
greater than or equal to these coefficients on a null data set, as above (Section 3.5). The 
correlation coefficient is weighted by uncertainties which are, for one data point, the square 
root of the sum of the squares of the uncertainties of the two quantities for which the 
correlation coefficient is being comput ed (see discussion on weighted least-squares fitting 



with uncertainties in both x and y by iBevingtonI Il969l ). Note that scaling all dust mass 



weights by the same scalar less than unity would not change the correlation coefficient 
between two dust mass fractions. The uncertainty weighting of the corr elation coe f ficient 



shows up in both numerator and denominator of the correlation coefficient IBevingtonI (119691 ) 
to the same power, so the scalar applied in the denominator would divide the scalar applied 
in the numerator to give 1 times the original correlation coefficient. Note that if r and the 
number, N, of data points that are being tested for correlation do not change, P will not 
change either because P depends only upon r and N. Values of r that significantly deviate 
from zero and values of P close to zero suggest significant correlation or anticorrelation. 
We deem as significant P<2.0%, resulting in |r| > 0.29. Stellar masses, disk-to-star mass 
ratios, ug-is and ni3_3i, and their uncertainties, plus multiplicity count s, stellar luminosities , 



and mass accretion rates used in this correlation anal ysis come from IWatson et al.l (120081 ). 



Submillimeter slopes were obtained from the study by Andrews fc Williamsl ( 20051 ). 



First we explore trends between dust mass fractions. All trends between crystalline 
components are positive. Figure 16 shows a positive correlation between warm enstatite and 
warm forsterite with correlation coefficient r=0.62 and probability P of having been drawn 
from a random distribution of 0.0%. This trend is fairly clear, with a heavy concentration 
of points at low mass fractions and a collection of points with higher enstatite and forsterite 
mass fractions. Warm enstatite also has a correlation with warm silica, with r=0.31 and 
P=1.2%. Figure 17 shows a positive correlation between cool forsterite and cool silica, with 
r=0.29 and P=2.0%. Cool forsterite also correlates with warm large grain mass fraction 
(r=0.32 and P=0.9%) and warm crystalline grain mass fraction (r=0.32 and P=1.0%). The 
correlation between warm silica and cool crystalline (cool enstatite plus cool forsterite plus 
cool silica) mass fractions is shown in Figure 18, with r=0.32 and P=0.9%. This correla- 
tion is likely related to the correlation between warm silica and cool forsterite (r=0.36 and 
P=0.4%). Warm enstatite also correlates with cool forsterite ( r=0.29 and P=1.9%). Coo l 



forsterite mass fraction and submillimeter slope {a reported by [Andrews fc Williamsl l2005l ) 
correlate, with r=-0.40 and P=2.0%. Lower a means flatter submillimeter slope (flux density 
at submillimeter wavelengths, Fiy^smm, is proportional to u"), implying that more growth of 
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grains to millimeter sizes (instead of sizes of a few microns) accompanies a greater abundance 
of cool forsterite. 

The general sense of these trends is that the crystalline dust abundances all correlate 
positively with each other. If the abundance of one type of crystal is high, it is likely that 
the abundances of the other two types will also be high. The gen eral finding that crysta l 
abundances track other crystal abundances has been noted before. iBouwman et al.l (l200ll ). 
via detailed modeling of the 10 fim complexes of Herbig Ae/Be stars, noted corr e lation 
between forsterite and silica abundance in Herbig Ae/Be stars, and Ivan Boekel et al.l (l2005l ) 
(using a similar analysis) noted correlation between enstatite abundance an d total crystalline 
abundance in a similar population of Herbig Ae/Be stars. IWatson et al.l (120081 ) measured 
crystalline abundances using indices computed from ratios of flux integrated over small 
wavelength bands characteristic of emission from pyroxene, olivine, and silica at 10 and 
33 fim (they call these bands Pio, Oio, Sio, and O33, respectively); in their study, they 
found warm pyroxene correlated with warm silica and warm olivine, and that warm olivine 



correlated with warm silica and cool olivine. From this we conclude, as did Watson et al. 



(l2008l ). that whatever produces the different kinds of crystals (forsterite, enstatite, and 
silica) produces them at a rate faster than the crystals can transform from one to another. 
As a specific example: if some silica forms by incongruent melting of enstatite, that amount 
of silica formed must be much less abundant than an y silica formed by another means 
(e.g., via incongruent melting of amorphous pyroxene; ISargent et al.l 120081 ). There is no 
significant correlation of crystalline silicate masses with disk-to-star mass ratio. A positive 
correlation of dis k-to-star mass ratio m ight be expected for crystalline dust produced by 
shock annealing (IHarker fc DeschI |2002| ). The lack of this correlation does not necessarily 
mean shock annealing is not connected to dust processing, however, as disk mass estimates 
are based on assumed opacities of dust at submillimeter wavelengths. The amount by which 
these assumed opacities differ from actual o pacities varies accord ing to the extent of grain 
growth to milhmeter sizes (see discussion by lD'Alessio et al.ll2006l ). 



Furlan et al.l (120061 ) computed and interpreted the SED continuum indices ng-ia and 



1113-25 as primarily measuring the degree of flaring of disks. Indices that are increasingly 
positive indicate increasingly flared disks, and more negative indices indicate increasingly 
flatter disks. The flattening of these disks was interpreted to mean that dust had s ettled 
from high in the disk atmosphere towards the disk midplane. IWatson et al.l (120081 ) used 
ni3_3i instead of ni3_25 to reduce contamination of the indices by the 20 /im amorphous 
silicate feature. We use the nfi_i3 and ni3_3i indices , interpreting them a s primarily measur- 
ing the degree of flaring of disks (IFurlan et al.ll2006l : IWatson et al.ll2008l ). We omit from our 
searches of correlation of parameters with dis k continuum indices t he points for the transi- 
tional and pre-transitional disks CoKu Tau/4 (iD'Alessio et al.ll2005l ). DM Tau and GM Aur 
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( ICalvet et al.ll2005l ). and LkCa 15 (lEspaillat et al.l 120071 ). These four objects have spectra 
whose continuum colors are very red due not to highly flared disks but, rather, to the clearing 
of almost all small dust grains in the innermost regions of the disks. The disk indices for 
these four systems would measure how clear they are of such dust in their inner regions, so 
they are not included in the search for trends of parameters with disk indices. 

Figure 19 illustrates a tight anticorrelation between warm crystalline grain abundance 
and Ug-is, with r=-0.50 and P=0.0%. This can be explained by correlations between each 
of the three warm crystalline types of grains and ug-ia. Warm forsterite and ug-is are 
anticorrelated with r=-0.42 and P=0.1%, warm enstatite and ng-is are correlated with r=- 
0.41 and P=0.1%, and warm silica and n(6 — 13) are correlated with r=-0.39 and P=0.2%. 
Cool enstatite correlates with ni3„3i (r=-0.31 and P=1.9%). 



These findings are similar to those of IWatson et al.l (120081 ) that greater warm olivine 
and warm silica, measured by Oio and Sio, respectively, anticorrelate with ng-ia and that 
greater warm enstatite, warm olivine, warm silica, and cool olivine (warm enstatite and 
cool olivine being measured by Pio and O33, respectively) anticorrelate with ni3_3i. The 
spectra described in Section 1.6 (low 10 /im complex equivalent width) all indicate, with 
the exceptions of DG Tau and V807 Tau, highly settled systems with abundant crystalline 
silicates, though with relatively large unc ertainties on these cryst alline abundances. The 
settling of disks is interpreted here, as did iFurlan et al. J2006h and IWatson et al.l (120081 ) , as 
part of the evolution of protoplanetary disks. The processing of dust in to crys t alline silicates 
is also interpreted as evolution of such disks. Further, the models by ICieslal (120071 ) predict 
crystalline mass fraction should be correlated with dust sedimentation both in inner and 
outer disk regions. Therefore, the correlation of crystalline abundance with disk settling is 
expected as a general result of protoplanetary disk evolution. 

Warm and cool large grain abundances both anticorrelate with ni3_3i. Figure 20 shows 
the correlation between warm large grain abundance and ni3_3i (r=-0.36 and P=0.5%). This 
trend is interpreted as an increase in large grain abundance expected as a result of evolution 
of protoplanetary disks. Figure 21 s hows that coo l large grain abundance also anticorrelates 
with ni3_3i (r=-0.32 and P=1.5%). iRettig et al.l (120061 ) found other tentative evidence for 
dust settling with grain growth. Grain growth is, in fact, expect ed to be the caus e of set thng 
of dust in protoplanetary disks. From first-principles modeling, IWeidenschillingI (119971 ) find 
in his simulations of grain growth and settling towards disk midplane of protoplanetary disks 
that at 60,000 years of evolution in the disk (the initial diameter of their dust grains was 1 
/xm), at 30 AU from the central star they find particles with diameters from one /im to near 
one centimeter. The range of ages assumed for Class II YSOs in the Taurus-Auriga star- 
forming region is 1-2 Myr, so substantial grain growth should have occurred for the 65 objects 
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whose spectra are analyzed in this study. Also, at a disk age of 0.1 Myr, IWeidenschilling 
( I1997I ) predicts that at a disk radius of 30 AU, the dust-to-gas mass ratio at 6 AU above 
the midplane should be 10~^ what it was at time zero, when the disk was initially well- 
mixed, at the same locati on in the disk. This rat io of the dust-to-gas mass ratio to its initial 
we ll-mixed valu e is w hat iD'Alessio et al.l (120061 ) call the settling parameter, e. As shown 
by iFurlan et al.l (120061 ) . the disk indices of some of the Class II YSOs in th e Taurus- Auriga 



star- fo rming region are consistent with such low values of e (0.01-0.001). iBouwman et al. 



(I2OO8I ) also find that the differing morphology of IRS spectra of YSOs are consistent with 
dust growth and settli ng. Large grain abund ance anticorrelating with ni3_3i is consistent 
with the expectation by IWeidenschillingI (119971 ) of increasing grain growth over time and also 
greater settling of larger dust grains towards the disk midplane. 

Also significant is the result shown in Figure 15a that more large gr ains are found in 



inner disk regions than outer disk regions. The observations and models by Ivan Boekel et al. 



(120041 ) required more large grains for disk regions inside 2 AU than fo r disk regions outside 
2AU f or spectra of two of the three Herbig Ae/Be stars they studied. [Andrews fc Williams 
(I2OO7I ) find that the surface density of protoplanetary disks decreases with radius in the disks 
as r~°'^ or r~^, so inner disk regions should be denser. This could encourage more rapid grain 
growth. 



However, as noted previously, the correlations between large grain abundances and 
ni3_3i are not very strong. One possible explanation for the weakness of correlation we find 
between warm large grain abundance and ni3_3i is that we are not measuring the true large 
grain abundance from each spectrum because we do not account for differing amorphous- 
to- crystalline component ratios in the large inhomogeneous aggregate grains from system to 
system. The more crystalline components the average aggregate large grain in a given disk 
has, the less reliable a measure our large amorphous silicate grains are of the large grain 
extent in such disks ( jMin et al.ll2008l ). 



Another possibility to explain the weakness of trend of warm large grain abundance with 
ni3_3i is that sufficiently high turbulen ce prevents rapid grain grow t h and settling of few- 
micro n-sized grains to the disk midplane (jPuUemond &: Dominikll2004j : iHubbard fc Blackman 
20061 ). Perhaps different disks can have the same degr ee of flattening but differing abun- 



dances of large grains in the optically thin surfa ce laye r. 



that collisions replenish small grains. However, ICieslal (120071 ) finds that turbulence can also 



Du 



lemond fc DominikI (120051 ) note 



encourage collisions that facilitate grain growth. These large grains, once formed, can be 
stored in a dead zone underneath the turbulent live zone, an d can be protected even if the 
dead zone is disturbed by the live zone above it (ICieslal 120071 ). 



Warm large grain abundance also correlates with the known number of stars in the star 
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system, with r=0.33 and P=0.7%, so perhaps multiphcity contributes to the dispersion seen 
in the trend between warm large grain abundance and ni3_3i. Figure 22 shows a comparison 
of the histogram of warm large grain abundance for single star systems and the same for 
multiple star systems - one immediately notices the higher large grain mass fractions for 
multiple systems than single systems. 

Finally, Figure 23 shows a fairly weak anticorrelation between warm large grain mass 
fraction and stellar mass, with r=-0.36 and P=0.5%. This may be related to the correlation 
between warm small amorphous grain mass fraction and stellar mass (r=0.40 and P=0.1%). 
In our sample, the spectra of low m ass stars indicate higher levels of warm large grain 



abundance than for higher-mass stars. lApai et al.l (120051 ) found in their study of six brown 



dwarf spectra higher mass fractions of crystalline grains for lower mass stars. They concluded 
that the region giving rise to the 10 /im feature was at smaller disk radii for less luminous, 
lower mass stars, and that in these innermost disk radii the crystalline dust abundance was 
higher. We interpret the trend of increasing warm (inner disk) large grain mass fraction with 
decreasing stellar mass similarly, in view of our finding of greater grain growth in inner disk 
regions than outer disk regions (Figure 15). 



6. Summary and Conclusions 

We have analyzed the dust composition of 65 T Tauri stars using spectra from the IRS 
on the Spitzer Space Telescope. These spectra show very prominent emission features and 
complexes at 10, 16, 19, 23, 28, and 33 fim wavelengths, which are characteristic of silica 
and silicate grains (both amorphous and crystalline) with sizes from submicron to a few 
microns. We have constructed spectral models for each of the spectra that include, firstly, 
blackbodies at two temperatures and, secondly, Planck functions at those two temperatures 
multiplied by scaled dust emissivities. These represent inner and outer disk emission from, 
firstly, the optically thick disk midplane and blackbody grains and, secondly, dust in the 
optically thin disk atmosphere with strong infrared resonances. The best fit for a given set 
of dust opacities is given by the pair of temperatures for which the global P^i' degree of 
freedom is a minimum. We have tested the opacities on high-quality IRS spectra to find 
the best opacities to represent emission from submicron grains of enstatite, forsterite, silica, 
and amorphous silicates and large (few micron size) grains in the protoplanetary disks of the 
Taurus-Auriga association. 

We conclude the following: 



High S/N spectra suggest the best crystalline silicate opacities to use in modeling 
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are those of iron-poor pyroxene (enstatite) and iron-poor olivine (forsterite). The 
best amorphous sihcate opacities to use are those whose iro n-to-magnesiuni rati os are 



nearly 1 (cosmic). This is consistent with the finding by iHarker et al.l (|2002| ) that 
cometary amorphous silicates require significant iron in their compositions in order to 
heat to inferre d dust temperatures . The adopted amorphous silicate closely matches 



the heating of iDraine fc Led (119841 ) "astronomical silicates". 



Though we use opacities of grains with 50% porosity and 5 fim radius of both amor- 
phous pyroxene and amorphous olivine, calculations of the covariance between these 
components in models of various high S/N spectra show they are highly degenerate. 
This means our models cannot readily distinguish between the two species as large 
grains, but our models are more sensitive to the sum of their masses. Submicron amor- 
phous pyroxene and submicron amorphous olivine are similarly degenerate but slightly 
less so. Cool silica is degenerate with cool submicron amorphous silicates as they share 
20 /im features. Cool enstatite and cool forsterite are somewhat degenerate because 
they share prominent 23, 28, and 33 /im features (only somewhat because the 28 fim 
feature of enstatite is relatively stronger than that of forsterite). 

A few spectra of T Tauri stars show evidence for extensive grain growth in that pro- 
toplanetary disks as indicated especially by their 10 fim complexes. 

Though most spectra are fit satisfactorily using our truncated CDE forsterite opacity, a 
few are not. The deviations of the 10, 23, and 33 /im features between our models and 
the spectra suggest olivine grains of greater size, greater porosity, greater iron content 
(namely, a nonzero iron content), or greater weighting towards extremely elongated or 
fiattened grain shapes in the shape distribution are required, or some combination of 
these. 

Most spectra that require silica are fit well by the annealed silica opacity in the standard 
dust model mixture. The spectra of FZ Tau and DN Tau suggest a silica poly morph 



other than annealed silica (cristobalite and tridymite; see iFabian et al.l l2000l ) would 
provide a better fit to their spectra. This would imply different cooling rates for silica, 
once f ormed, than implied by the presence of the higher temperature polymorphs of 



silica (ISargent et al.ll2008l ). However, the 12.6 /im features of these two exceptional 
spectra are only mildly inconsistent with the annealed silica in the dust models, so the 
suggestion of an alternative silica polymorph is weak. 

The models of the best enstatite exemplars, FN Tau and DH Tau, do not fit very well 
the strong crystalline silicate emission peaks at 28 and especially 33 /im, suggesting 
enstatite or forsterite dust at temperatures significantly lower than those used by our 
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models is required by the data. It is unknown why such cold material is required to fit 
spectra indicating high enstatite abundance. 

Transitional and pre-transitional disks in this sample require negligible crystalline sil- 
icates and modest amounts of large grains (DM Tau being the exception, requiring 
significant large grains) to fit their spectra relative to the rest of the disks whose spec- 
tra were analyzed in this study. A similar (small) number of spectra of systems not 
known to be transitional or pre-transitional also require high abundances of submicron 
amorphous silicates. 

A few spectra have 10 fim complexes that have small equivalent widths. Most of these 
are fit well by crystalline silicates as opposed to large grains, which is inconsistent with 
the idea that decreasing peak-to-continuum ratio of the 10 /im complex necessarily 
implies grain growth. 

The goodness of fit of our models to the majority of our spectra means the opacities of 
large crystalline silicate grains are not typically required. This, in turn, implies that, 
if grains grow to larger sizes by agglomeration of grains of differing composition, then 
on average in the Taurus-Auriga YSO population the more abundant subcomponents 
of these heterogeneous aggregate large grains are amorphous rather than crystalline. 

Higher abundances of large grains are found in inner disk regions than outer disk 
regions, suggesting grain growth occurs more rapidly in inner disk regions. 

Crystalline silicate abundances are very similar in the inner and outer disk regions. 



This contrasts interestingly with the finding by Ivan Boekel et al.l (120041 ) that regions 
inside 2 AU in disks around Herbig Ae/Be stars have much more processed dust than 
regions outside of 2 AU. F04147+2822 exemplifies the contrast, having prominent 19, 
23, 28, and 33 fim complexes but very little indication of crystallinity based on its 
10 fim complex. This suggests the inner region of this disk giving rise to the 10 /im 
complex has a lower abundance of crystals than the outer region giving rise to emission 
past 19 fim wavelength. 

Cool forsterite also correlates with flatter submillimeter SEDs, suggesting grain growth 
to millimeter sizes is more extensive for disks with more cool forsterite. 

Each crystalline silicate abundance always correlates positively and significantly with 
abundances of other crystalline silicate species. This suggests that whatever produces 
crystalline silicates does so at a faster rate than that at which any of the crystalline 
silicate species transform into any of the others. This, in turn, suggests that amorphous 
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dust is processed into crystalline dust at a rate greater than one kind of crystalline 
silicate species can transform into another. 

• Crystalline silicates in inner disk regions are more abundant for bluer, flatter disks, 
those with more advanced settling of dust grains towards disk midplane. 

• Large grains in the inner disk regions are, on average, more abundant for more settled 
(bluer) disks. There is a fair degree of dispersion in the inner disk large grain abundance 
for a given degree of settling, suggesting different rates of settling of large grains in 
different disks. There is also an indication that the abundance of large grains in the 
inner disk correlates with known multiplicity of stellar system. 

• The spectra of disks around less massive stars indicate higher abundances of large 
grains in their inner disk regions. 
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A. Error Bars and Mispointing 



We corrected for slight mispointing of t he telescope from th e standard nod positions by 
scahng the mispointed nods, as described by lSargent et al.l (120081 ). For observations taken by 
the low spectral resolution modules SL and LL, only mispointing in the dispersion direction 
(perpendicular to the length of the slit) could give rise to loss of flux. Mispointing of the high 
resolution slits along the length of the slit (spatial direction) towards the center of the slit 
can actually result in obtaining more flux than is present if one does not consider the RSRFs 
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used to calibrate flux were constructed from observations of flux standards at nominal nod 
positions. More flux enters through the slit if the target is positioned closer to the center of 
the slit than a nominal nod position (like the flux calibrator is); in this case, one must apply 
a scalar factor less than 1 to the affected nod spectrum to correct for the extra fraction of 
total signal from the target not measured for the calibrator. SH has a smaller slit and is 
therefore more sensitive to mispointing in either spectral or spatial directions; therefore, the 
relative amount of scalar correction is typically higher for SH than LH. The scalars used are 
given in Table Al. 



B. Test of Large Grains with Heterogeneous Composition 



We conducted tests of the combination of Bruggeman EMT and Mie Theory used to 
compute effec t ive co mplex dielectric constants and opacities for our large grains. Note that 
Sargent et al.l (120061 ) found that opacity curves of sub-micron spherical 60% vacuum (porous) 
grains of forsterite and a-quartz were identical to those computed for solid grains of the 
same materials in the CDE shape distribution. We have computed the effective complex 
dielectric function for material that is, by volume, 60% vacuum (c omplex dielectric fu nction 



equal to unity), 20% forsterite (complex dielectric function by ISogawa et al 



2006 land 



20% amorphous olivine MgFeSi04 (complex dielectric function by lDorschner et al.lll995l ) as 
follows. First we computed using Bruggeman EMT the effective complex dielectric function 
of 50% amorphous olivine and 50% forsterite a-axis; we also computed the same effective 
complex dielectric functions for 50/50 mixtures of amorphous olivine and each of the other 
two forsterite crystallographic axes, b and c. Then we computed using Bruggeman EMT 
the effective complex dielectric constants for material that is 60% by volume vacuum and 
40% the 50/50 mixture of amorphous olivine and forsterite; again, the same was done for 
the effective complex dielectric functions involving forsterite b and c axes. This resulted in 3 
sets of complex dielectric functions, one for each of the 3 crystallographic axes of forsterite, 
of material that is, by volume, 60% vacuum, 20% forsterite, and 20% amorphous olivine. 
Next we computed opacities for each of the three sets of effective complex dielectric functions 
for Rayleigh-limit-size spheres of these materials. Finally, we averaged the three opac ities 



together to account for anisotropy of optical constants (see iBohren fc Huffman! Il983l ) for 
forsterite. We then calculated the flux density from 1 lunar mass of such grains 140 parsecs 
away at 300K. 1% relative uncertainties were assumed, which are about the best we can 
expect from our sample of 65 spectra. 



The 7.7-37 /im part of this spectrum was then modeled using our standard opacities, 
which are listed in the notes for Table 6, with the exception of using the opacities of forsterite 
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and amorphous olivine grains in the CDE shape distribution (instead of tCDE and CDE2). 
We used the CDE for these two d ust types in anticipation of the porous-CDE correspondence 



reported by ISargent et aLl (120061 ) . In Figure Al, we show the fit of the model to this test 
spectrum. Although none of the dust opacities were precluded at the outset, the only 
opacities used were the CDE opacities for forsterite and amorphous olivine, with masses 
of 0.45 and 0.43 Mmoou, respectively. The fit is very good to all wavelengths, though we 
note deviations for the strongest features at 11.4, 19, and 23 /im; we also note that the 



correspondence between porous grain and CDE opacities reported by I Sargent et al.l (120061 ) 
was not perfect. The model used grains at 308K, only slightly different from the assumed 
300K. This slightly higher than assumed temperature explains the slightly lower (0.88 Mmooti) 
than assumed (1 Mmooti) total mass of dust. We take this as support for the conclusion by 



Min et al.l (120081 ) that opacities of aggregates can be expressed as a sum of the opacities 
of homogeneous grains composed of the same materials as the aggregate's components, as 
applied to heterogeneous aggregates in the Rayleigh size limit. 

Next we conducted a similar test, but for large aggregates. We used the same effec- 
tive complex dielectric function of vacuum, forsterite, and amorphous olivine in a 60/20/20 
volume ratio as before but this time to compute opacities for 5 /xm radius grains of the 
material. We computed the test spectrum assuming 1 lunar mass of grains at 300K located 
140pc away. Here, we used our standard opacities exactly as hsted in the notes for Table 6, to 
gauge the effect of large grains in protoplanetary disks being aggregates of both amorphous 
and crystalline grains. Dust at only one temperature was allowed, but no dust species at 
this temperature was precluded from the model at the outset. This time, the best-fit model 
dust temperature is lower, 264K, and enstatite is used in addition to forsterite and large 
(5 fim radius porous) grains of amorphous olivine, though at a lower significance level than 
either the amorphous olivine or forsterite. We show this model in Figure A2. The total 
mass of the dust used by the model, 1.25 lunar masses, is a sum of 0.96 Mmoou of large 
amorphous olivine plus 0.20 Mmooti of forsterite, 0.09 Mmoou of enstatite, and 0.001 Mmooti 
of submicron amorphous pyroxene, which is somewhat greater than the mass assumed for 
the aggregates. This increase in mass was needed to compensate for the lower tempera- 
ture. The model becomes an increasingly poorer fit at increasing wavelengths. Part of this 
failure may be due to the use of the restricted shape distribution (tCDE) for the forsterite 
instead of the more extreme shape distribution CDE (used previously to test the Rayleigh- 
limit aggregate). An additional population of forsterite grains at a lower temperature than 
used may improve the fit. This suggests any serious failures of our standard model to fit 
spectra with prominent crystalline silicate features may imply a significant amount of grain 
growth of aggregate grains with significant crystalline silicate abundances (here, 50% of the 
solid mass of the aggregate). This test also implies that the crystalline abundance could be 
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underestimated by 50% (0.29/1.25 instead of 0.625/1.25 fractional abundance of crystalline 
silicates in heterogeneous aggregates) of the true value, and it implies that the large grain 
abundance could be slightly underestimated (0.96 instead of 1.25 lunar masses of large (5 
/im radius) grains) by 100 x (0.29/1.25)% = 23%. 



C. Importance of Large Grain Opacity 

We tested the necessity of using opacities of large grains in our models. We used our 
standard dust model on IS Tau, one of the objects whose spectrum was modeled using 
a significant amount of large grains, but we eliminated the opacities of large amorphous 
pyroxene and large amorphous olivine in the fit. The resulting model, shown in Figure 
A3, is a poorer fit to the 10 /xm complex, not "filling out" the full width of the 10 fim 
complex. Reduced ^ot this large-grains-precluded model was 4.9, which is 1.6 higher than 
the reduced of 3.3 for the model with large grains. We conclude the opacities of large 
amorphous grains are required in our models. 



D. Porosity Test 



We tested how sensitive our models are to the exact combination of porosity and grain 
size when we compute and use opacities of large grains in our models to detect grain growth. 
Instead of 5 /im radius 60% porous grains, we computed opacities for 20 /im radius grains 
that are 88% porous. We chose this s i ze to preserve the optical depth through the center of 
the grain. As noted by lSargent et al.l (120061 ) . 5 fim radius 60% porous grains were chosen to 
preserve the optical depth throu gh the center of the 2 /im solid grains of amorphous olivine 
used by iBouwman et al.l (1200 ll ) to represent grain growth in modeling spectra of Herbig 
Ae/Be stars. This preservation of the optical depth through the center of the grain resulted 
in a close correspondence of the opacities of 5 fj,m radius porous grains and of 2 /im radius 
solid grains. Although 90% vacuum (by volume) grains would be the logical extension 88% 
vacuum 20 /im radius grains proved to result in an opacity most like that of 5 /im radius 
60% vacuum grains. As a test of this opacity, we modeled the spectrum of UZ Tau/e, noted 
in Section 4 to require a large amount of large grains; here, we replace the opacity of 5 /im 
60% porous grains of amorphous olivine with that of 20 /im 88% porous grains of amorphous 
olivine. We show the best-fit model in Figure A4, compared to the model shown in Figure 4 
and described in Section 4. The fit is not much different, and the reduced is 5.6 instead of 
the 5.4 obtained with our standard opacities. Therefore, we cannot make precise statements 
on grain size and porosity individually for a given protoplanetary disk based only on the IRS 
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spectra. 
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Fig. 1. — Model of the spectrum of the Trapezium. Model components are at bottom. 
Blue lines represent model components at higher temperature, red lines represent model 
components at lower temperature. Solid lines are blackbodies, dotted lines are submicron 
amorphous pyroxene, dash-double-dotted lines are submicron amorphous olivine, dot-long- 
dash lines are large amorphous pyroxene, dot-double-dash lines are large amorphous olivine, 
dot- short- dash lines are enstatite, long-dashed lines are forsterite, and short-dashed lines are 
silica. 
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Fig. 2. — Contours of per degree of freedom for models of the Trapezium spectrum with 
dust at the same temperatures as the model shown in Figure 1, but with dust weights allowed 
to be negative (the "negative-allowed" model). The first three contours of reduced are 
labeled. The "x" gives the best-fit values of cool small amorphous pyroxene and cool small 
amorphous ohvine for this model, resulting in reduced of 1.4. The "+" gives the best- fit 
values of the same two dust components for the standard model (the "negative-not-allowed" 
model) shown in Figure 1, resulting in reduced of 1.6. This value, 1.6, is lower than that 
indicated on the plot, 2.0, because in the negative-not-allowed model, the mass weights of 
the other dust components were allowed to vary to obtain the global minimum in reduced 
X^, while in the negative-allowed model they were not. Note that nearly the same minimal 
level of reduced can be obtained by replacing a given amount of one of the two dust 
components by about the same mass of the other. 
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Fig. 3. — Histograms of the correlation coefficients obtained for all 65 spectra of Taurus- 
Auriga T Tauri stars analyzed in this study, for four representative pairs of dust components. 
Varying degrees of degeneracy, or anticorrelation, of components are seen, with the most 
degeneracy coming from pairs of amorphous silicate grains at the same temperature. 
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X (|Lim) 

Fig. 4. — Spectra indicating substantial inner disk grain growth. Model components for 
model of UZ Tau/e are at bottom; blue lines represent model components at higher tem- 
perature, red lines represent model components at lower temperature. Solid lines are black- 
bodies, dotted hues are small (submicron) amorphous pyroxene, dash-double-dotted lines 
are small (submicron) amorphous olivine, dot-long-dash lines are large amorphous pyrox- 
ene, dot-double-dash lines are large amorphous olivine, dot-short-dash lines are enstatite, 
long-dashed lines are forsterite, and short-dashed lines are silica. 
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Fig. 5. — Spectra with prominent forsterite features. Model components are shown for IS 
Tau, same style and color convention as Figure 4. 
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Fig. 6. — Spectra indicating forsterite, but not fit well by our forsterite opacity. Model 
components are shown for FX Tau, same style and color convention as Figure 4. 
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Fig. 7. — Spectra of DN Tau, ZZ Tau, and FZ Tau. DN Tau and FZ Tau require silica but 
are not fit well by our silica opacity, while ZZ Tau requires silica and is fairly fit well by our 
silica opacity. Model components are shown for FZ Tau, same style and color convention as 
Figure 4. 
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Fig. 8. — Spectra with prominent enstatite features. Model components are shown for FN 
Tau, same style and color convention as Figure 4. 
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Fig. 9. — Spectra whose models require high abundances of submicron amorphous sihcates. 
Model components are shown for FM Tau, same style and color convention as Figure 4. 
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Fig. 10. — Spectra whose models require fairly high abundances of submicron amorphous 
silicates, though not as high as those in Figure 9. Model components are shown for HQ Tau, 
same style and color convention as Figure 4. 
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Fig. 11. — Spectra with low equivalent width 10 fim complexes. Model components are 
shown for DF Tau, same style and color convention as Figure 4. 
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Fig. 12. — Spectra whose models indicate mixed compositions, with no dust component 
clearly dominant over another. Model components are shown for FP Tau, same style and 
color convention as Figure 4. 
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Fig. 13. — Spectra whose models indicate mixed compositions, with no dust component 
clearly dominant over another. Model components are shown for GI Tau, same style and 
color convention as Figure 4. 
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Fig. 14. — Spectra whose models indicate mixed compositions, with no dust component 
clearly dominant over another. Model components are shown for VY Tau, same style and 
color convention as Figure 4. 
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Fig. 15. — Warm versus cool dust histograms. In the upper left is the large grain mass 
fraction histogram. In the middle left is the crystalline grain mass fraction histogram. In 
the lower left is the enstatite grain mass fraction histogram. In the upper right is the 
forsterite grain mass fraction histogram. In the lower right is the silica grain mass fraction 
histogram. The number of models requiring warm grain mass fractions within a given bin are 
represented by the height of the blue rectangle to the left within that bin. Orange rectangles 
to the right are for cool grain mass fractions. 
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Fig. 16. — Warm forsterite mass fraction versus warm enstatite mass fraction trend plot. 
The thick sohd hne running through the points is the trendhne consistent with the computed 
correlation coefficient, indicated on the plot. 
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Fig. 17. — Cool forsterite mass fraction versus cool silica mass fraction trend plot. The 
thick solid line running through the points is the trendline consistent with the computed 
correlation coefficient, indicated on the plot. 
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Fig. 18. — Cool crystalline mass fraction versus warm silica mass fraction trend plot. The 
thick solid line running through the points is the trendline consistent with the computed 
correlation coefficient, indicated on the plot. 
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Fig. 19. — Warm crystalline mass fraction versus ng-is trend plot. The thick solid line run- 
ning through the points is the trendline consistent with the computed correlation coefficient, 
indicated on the plot. 
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Fig. 20. — Cool crystalline mass fraction versus ni3_3i trend plot. The thick solid line running 
through the points is the trendline consistent with the computed correlation coefficient, 
indicated on the plot. 
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Fig. 21. — Warm large grain mass fraction versus ni3_3i trend plot. The thick solid line run- 
ning through the points is the trendline consistent with the computed correlation coefficient, 
indicated on the plot. 
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Fig. 22. — Warm large grain mass fraction histograms for single systems (blue) and multiple 
systems (green). 
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Fig. 23. — Warm large grain mass fraction versus stellar mass trend plot. The thick solid 
line running through the points is the trendline consistent with the computed correlation 
coefficient, indicated on the plot. 
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Fig. Al. — Model of Rayleigh-limit size porous aggregate grain composed of amorphous 
olivine and forsterite. Style and color convention of model components same as for Figure 4. 
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Fig. A2. — Model of 5 fim radius porous aggregate grain composed of amorphous olivine 
and forsterite. Style and color convention of model components same as for Figure 4. 
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Fig. A3. — Model of IS Tau, without using large grains at either model dust temperature. 
Reduced is 4.9, which is 1.6 higher than the 3.3 computed for the model of IS Tau using 
large grains (see Fig. 5). Style and color convention of model components same as for Figure 
4. 
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Fig. A4. — Model of UZ Tau/e, replacing 5 /xm radius 60% vacuum porous grains of amor- 
phous olivine with 20 /im radius 88% vacuum porous grains of amorphous olivine. Style and 
color convention of model components same as for Figure 4. 
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Table 1: Significance (mass/uncertainty) of Opacities for Forsterite Exemplars 



Opacity 


DK Tau 


GN Tau 


HBC 656 


IS Tau 


ROXs 42C 


V836 Tau 


V955 Tau 




CfoJ98 





4.3 


5.0 


1.4 


0.7 








< 


ColJ98 


4.5 






















ChoJ98 




















0.1 




CfaJ98 














0.9 










CFo91.4 

























CFo89.3 

















1.6 







CF086.5 




















1.8 




CFo84.4 


0.4 






















CF08O 

























CF02I.8 








0.6 
















Cannfor 

























Cfohyb 








1.8 


1.0 








0.3 




CS6CDE 








0.9 


0.1 








0.8 




CS6CDE2 


1.3 





2.6 


2.7 


3.2 


1.1 


2.2 




WfoJ98 

























WolJ98 

























Who.I98 














0.3 










WfaJ98 


3.1 





0.9 


2.5 


2.4 


0.4 


1.7 




WF09I.4 




















1.2 




WF089.3 

























WF086.5 











1.3 













WFo84.4 


0.1 






















WF08O 





1.2 


1.1 





0.6 










WF02I.8 








2.1 
















Wannfor 


4.0 





2.6 


1.2 





3.3 







Wfohyb 

















0.6 







WS6CDE 





0.1 


0.9 








1.2 







WS6CDE2 


0.6 


0.4 


2.6 


0.5 


0.9 





0.4 





Note. — A given entry in a column is the mass weight of the component named in the first entry on the 
same row divided by its uncertainty, giving the significance of a particular opacity used in a given model. The 
naming convention is as described in Section 3.2. Arrows in the final column indicate recurri ng opacities as 



discus sed in the text. The opacity used most at both temperatures is the CDE2 opacity using ISogawa et al 



(|2006j) optical constants, with the CDE opacity from the same optical constants also frequently included. 
Both of these are indicated by double-lined arrows. 
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Table 2: Significance (mass/uncertainty) of Opacities for Enstatite and Forsterite Exemplars, 
Jointly 



Opacity 


DH 
Tau 


FN 
Tau 


Haro 
6-37 


HK 

Tau 


DK 

Tau 


GN 
Tau 


HBC 
656 


IS 

Tau 


ROXs 
420 


V836 
Tau 


V955 
Tau 


Cen90 








1.6 





1.9 


2.1 











0.9 





CenSO 



































ConTO 




















1.0 














CcnSO 




















0.1 














CenOO 











0.4 











1.0 


2.1 





1.4 


CJ98CDE2 


2.7 


2.7 





0.7 





1.5 


2.5 


0.7 


0.9 








CKOso 



































CJ98b 














1.4 




















CJ98c 


0.1 

















3.4 





2.7 


0.4 





CS6CDE2 


























2.3 








CSetCDEO.lO 








0.9 


2.6 


2.2 








4.6 





2.4 


3.6 


CSetCDEO.OB 


























1.5 








CS6tCDE0.06 

















2.6 

















CS6tCDE0.04 


4.0 


3.3 














8.1 














Wcn90 


6.2 


3.5 


3.6 


3.6 


3.5 


0.2 


4.9 


2.1 


0.3 


2.5 


1.3 


WcnBO 


























0.3 





0.5 


WonTO 


0.7 
































WonSO 



































WonOO 


1.8 





2.6 


1.9 


2.4 





3.2 


1.0 





4.3 





WJ98CDE2 








0.3 


0.2 








0.5 














WKOso 





2.0 





























WJ98b 



































WJ98c 



































WS6CDE2 



































WSetCDEO.lO 


4.2 


0.9 


1.2 





2.9 





3.2 


1.4 





2.0 





WSetCDEO.OB 




















2.4 





1.3 








WS6tCDE0.06 
































1.2 


WS6tCDE0.04 











0.9 





1.3 


















Note. — Entries in the table are mass weight of a component divided by uncertainty, as in Table 1. Arrows 
in the final column indicate recurring opacities as discussed in the t ext. The enstatite opacity used most at 
both temperatures is the "en90" opacity from Chihara et al. ( 20021 ). in dicated by double-lin ed arrows. The 
optimal forsterite opacity is the tCDE opacity using optical constants bv lSogawa et al.l (|2006| ) with bounding 
parameter t=0.10, indicated by double-lined arrows. 
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Table 3: Significance (mass/uncertainty) of Opacities for Amorplious Silicate Exemplars 



Opacity 


CoKu Tau/4 


DM Tau 


FM Tau 


GM Aur 


LkCa 15 


TW Cha 


UY Aur 




CD9501Sm 























^= 


CD9501Lg 

























CD7901Sm 


5.5 






















CD7901Lg 

























CD95Py5Sm 


12.4 

















1.4 




CD95Py5Lg 

























CD95PylOSm 


6.6 





7.3 


14.1 











< 


CD95PylOLg 

























CJ94PySm 











4.2 


12.0 


10.6 







CJ94PyLg 

























CD79PySm 


0.8 


3.0 





7.5 





0.4 


5.4 


< 


CD79PyLg 

























CSD96PySm 

















7.0 







CSD96PyLg 

























WD9501Sm 


8.9 


5.8 


15.0 


16.2 


23.2 


22.4 


9.3 




WD9501Lg 

























WD7901Sm 





4.3 





0.4 


1.1 





1.7 




WD7901Lg 





1.0 











2.5 


0.1 




WD95Py5Sm 





3.0 



















WD95Py5Lg 





15.1 


0.3 
















WD95PylOSm 





0.5 


0.5 











2.0 




WD95PylOLg 








3.2 
















WJ94PySm 

























WJ94PyLg 














2.5 


4.7 







WD79PySm 

























WD79PyLg 





1.7 


1.1 











1.8 




WSD96PySm 











7.7 


3.5 





1.8 




WSD96PyLg 


























Note. — Entries in the table are mass weight of a component divided by uncertainty, as in Table 1. Arrows 
in the final column indicate recurring opacities as discussed in the text. Of the am orphous olivines, t he on e 
included in the most models at the highest significance is the amorphous olivine from lPorschner et al. (1995), 
indicated by double-lined arrows. Of the amorphous pyroxenes, the one that s hows up in the mos t models at 
the highest significance is the amorphous pyroxene of cosmic composition bv I Jaeger et al.l (|1994| ). indicated 
by double-lined arrows. 
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Table 4: Trapezium Dust Model Parameters 



Dust 


281 K 






145 K 






Species 


Value 


Unc. 


Sgnfcnc. 


Value 


Unc. 


Sgnfcnc. 


(1) 


(2) 


(3) 


(1) 


(5) 


(6) 


(7) 


Blackbody 


59.6 


8.8 


6.8 








Sm. Am. Py. 


1.90 


0.63 


3.0 








Sm. Am. 01. 


2.05 


0.75 


2.7 


358 


47 


7.6 


Silica 


0.367 


0.292 


1.3 


23.0 


29.9 


0.8 



Note. — Col. (1): Dust component. "Sm. Am. Py." and "Sm Am. 01." denote small amorphous 
pyroxene and small amorphous olivine, respectively. Cols. (2) and (5): Parameter value for warm and cool 
dust, respectively. Cols. (3) and (6): Uncertainties on parameter values for warm and cool dust, respectively. 
Cols. (4) and (7): Significances of component, being the parameter value divided by its uncertainty, for warm 
and cool dust, respectively. Parameter values and uncertainties for blackbody are in units of 10~^^ steradians, 
while parameter values and uncertainties for other dust components are in units of 10~^^ indicates 
the component at the temperature indicated by placement in the columns was not included in the model, so 
its mass is zero. Other components not included in the model (also having zero mass) are large amorphous 
pyroxene, large amorphous olivine, enstatite, and forsterite. per degree of freedom for this model is 1.6. 



Table 5: Degeneracies between Dust Components 



Component 



Component 



Corral. 
Coeff. 



(1) 




(3) 



Wm. Sm. Am. Py. 

Wm. Lg. Am. Py. 

Wm. Enstatite 



Wm. Sm. Am. 01. 
Wm. Lg. Am. 01. 



Wm. Silica 



-0.73 
-0.72 

-0.50 



CI. Sm. Am. Py 
CI. Lg. Am. Py 
CI. Enstatite 
CI. Enstatite 



CI. Sm. Am. 01. 
CI. Lg. Am. 01. 

CI. Forsterite 



CI. Silica 



0.82>r>-0.89 
0.87>r>-0.89 

-0.46 
0.29>r>-0.50 



Note. — Cols. (1) and (2): pair of dust components with significant degeneracy. "Wm." denotes warm, 
"CI." denotes cool, "Sm." denotes small, "Lg." denotes large, "Am." denotes amorphous, "Py." denotes 
pyroxene, and "01." denotes olivine. Col. (3): Correlation coefRcients between pair of dust components 
obtained from covariance matrices. The probability, P, of the correlation coefBcient between each of the 
seven pairs listed in the table being drawn from a random distribution of data is less than or equal to 
0.1%, assuming there are 133 independent data points in each spectrum for which the a covariance matrix 
is computed. 



Table 6. Dust Model Parameters 









small 


small 


large 


large 








Total 






Temp 




Amorphous 


Amorph. 


Amorph. 


Amorph. 


Crystalline 


Crystalline 


Crystalline 


Dust 




Object 


(K) 




Pyroxene'' 


Olivine'' 


Pyroxene'^ 


Olivine'' 


Enstatite'' 


Forsterite' 


Silicas 


Mass 


2 

X 

d.o.f. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


AA Tau 


127 


123.8 ± 3.9 


75.1 ± 30.1 


0.0 ± 18.8 


0.0 ± 21.6 


0.0 ± 14.3 


6.2 ± 11.7 


7.6 ± 10.6 


11.1 ± 14.4 


59.1 


2.9 




545 


1.22 ± 0.02 


9.7 ± 4.4 


28.9 ± 5.0 


34.2 ± 6.1 


22.6 ± 5.8 


4.5 ± 3.5 


0.0 ± 3.2 


0.0 ± 2.5 


0.937 




BP Tau 


103 


64.3 ± 10.9 


15.5 ± 4.9 


24.4 ± 3.8 


45.2 ± 4.6 


0.0 ± 2.6 


8.4 ± 2.5 


3.0 ± 2.1 


3.5 ± 3.7 


981 


3.6 




431 


2.73 ± 0.05 


45.6 ± 6.4 


27.7 ± 5.2 


9.3 ± 6.9 


0.0 ± 6.4 


12.9 ± 3.8 


4.6 ± 3.4 


0.0 ± 2.8 


1.95 




CI Tau 


115 


312.4 ± 9.93272 


80.2 ± 17.9 


11.8 ± 10.0 


0.0 ± 11.6 


0.0 ± 7.6 


0.0 ± 6.9 


8.0 ± 6.3 


0.0 ± 9.2 


278 


3.3 




545 


1.50 ± 0.03 


0.0 ± 3.5 


35.1 ± 4.1 


47.6 ± 4.9 


11.4 ± 4.6 


4.7 ± 2.7 


0.7 ± 2.5 


0.5 ± 1.9 


1.66 




CoKu Tau/3 


103 


0.0 ± 0.0 


0.0 ± 7.8 


3.0 ± 5.6 


24.9 ± 5.5 


54.0 ± 7.6 


3.9 ± 3.6 


9.1 ± 3.4 


5.2 ± 5.4 


488 


4.7 




602 


0.975 ± 0.019 


11.1 ± 3.2 


13.2 ± 3.4 


28.0 ± 4.4 


36.6 ± 4.3 


7.0 ± 2.6 


4.1 ± 2.4 


0.0 ± 1.8 


1.06 




CoKu Tau/4 


114 


354.2 ± 8.2 


32.8 ± 3.9 


47.6 ± 4.2 


13.8 ± 3.7 


0.0 ± 2.8 


0.0 ± 2.3 


2.3 ± 2.0 


3.4 ± 2.7 


677 


2.2 




200 


0.420 ± 0.230 


8.4 ± 10.0 


91.6 ± 22.4 


0.0 ± 14.1 


0.0 ± 12.5 


0.0 ± 6.9 


0.0 ± 5.1 


0.0 ± 7.6 


4.64 




CW Tau 


151 


142.2 ± 3.9 


0.0 ± 26.8 


24.8 ± 17.8 


0.0 ± 21.1 


0.0 ± 14.0 


15.5 ± 12.7 


29.8 ± 15.7 


29.9 ± 17.3 


59.1 


5.3 




659 


1.52 ± 0.03 


0.0 ± 7.5 


51.6 ± 11.7 


42.1 ± 10.8 


0.0 ± 12.4 


2.9 ± 5.9 


0.0 ± 5.2 


3.4 ± 4.2 


0.643 




CX Tau 


151 


58.3 ± 1.4 


27.2 ± 15.1 


62.7 ± 19.8 


0.0 ± 16.0 


0.0 ± 11.5 


0.0 ± 8.8 


3.4 ± 7.5 


6.8 ± 9.4 


27.0 


2.6 




602 


0.400 ± 0.013 


7.3 ± 6.8 


15.2 ± 7.4 


24.4 ± 9.9 


38.9 ± 9.9 


7.6 ± 5.4 


3.2 ± 5.1 


3.4 ± 4.2 


0.268 




CY Tau 


127 


22.4 ± 1.0 


22.6 ± 37.7 


0.0 ± 30.1 


0.0 ± 34.0 


0.0 ± 21.9 


0.0 ± 20.2 


48.2 ± 38.1 


29.2 ± 28.6 


9.50 


4.0 




602 


0.319 ± 0.006 


0.0 ± 22.1 


69.9 ± 39.9 


0.0 ± 33.1 


0.0 ± 30.1 


12.8 ± 17.0 


0.0 ± 13.839 


17.3 ± 15.1 


0.0460 




DD Tau 


127 


220.2 ± 8.4 


80.5 ± 18.1 


0.0 ± 10.7 


0.0 ± 11.8 


0.0 ± 7.9 


10.1 ± 6.8 


9.3 ± 6.4 


0.0 ± 9.0 


233 


3.5 




488 


4.72 ± 0.07 


26.3 ± 9.1 


0.0 ± 9.7 


7.9 ± 13.0 


46.2 ± 12.3 


13.7 ± 6.8 


5.8 ± 6.1 


0.0 ± 4.8 


1.42 




DE Tau 


115 


226.5 ± 7.7 


84.7 ± 9.9 


5.2 ± 5.6 


0.0 ± 6.2 


0.0 ± 4.0 


0.0 ± 3.5 


10.1 ± 3.4 


0.0 ± 4.6 


390 


4.7 




545 


1.17 ± 0.02 


0.0 ± 6.1 


59.2 ± 10.4 


0.0 ± 10.1 


28.0 ± 8.2 


5.0 ± 4.8 


7.8 ± 4.3 


0.0 ± 3.4 


0.600 




DF Tau 


163 


46.6 ± 2.2 


0.0 ± 14.1 


0.0 ± 10.2 


0.0 ± 12.3 


86.4 ± 21.3 


0.0 ± 6.8 


4.6 ± 5.8 


9.0 ± 8.0 


56.2 


5.3 




716 


1.34 ± 0.02 


63.6 ± 27.0 


0.0 ± 17.4 


0.0 ± 23.8 


0.0 ± 23.4 


18.5 ± 12.5 


17.9 ± 11.6 


0.0 ± 7.6 


0.250 




DG Tau 


127 


2810.9 ± 56.2 


86.3 ± 18.4 


0.0 ± 11.0 


0.0 ± 11.3 


0.0 ± 7.3 


0.0 ± 5.8 


11.7 ± 6.4 


1.9 ± 8.2 


1350 


6.4 




431 


23.1 ± 0.4 


0.0 ± 335.1 


0.0 ± 531.6 


0.0 ± 745.9 


0.0 ± 855.8 


0.0 ± 326.0 


100.0 ± 1347.6 


0.0 ± 162.1 


0.113 




DH Tau 


127 


89.4 ± 2.8 


6.7 ± 14.5 


48.9 ± 13.4 


0.0 ± 12.4 


0.0 ± 8.5 


17.9 ± 7.7 


26.6 ± 8.9 


0.0 ± 7.0 


68.3 


8.4 




944 


0.0953 ± 0.0032 


0.0 ± 4.1 


7.5 ± 4.9 


0.0 ± 7.4 


47.3 ± 6.7 


26.0 ± 4.5 


14.8 ± 3.8 


4.3 ± 2.1 


0.115 




DK Tau 


175 


74.0 ± 2.1 


0.0 ± 14.0 


9.1 ± 9.4 


19.0 ± 11.3 


15.9 ± 8.5 


17.1 ± 7.4 


15.4 ± 6.8 


23.4 ± 9.1 


54.3 


5.7 




830 


0.631 ± 0.020 


0.0 ± 2.2 


1.5 ± 2.6 


30.3 ± 3.1 


54.2 ± 3.6 


7.1 ± 1.9 


4.8 ± 1.8 


2.0 ± 1.1 


1.58 




DL Tau 


139 


147.2 ± 5.2 


0.0 ± 25.1 


0.0 ± 17.3 


36.0 ± 20.3 


41.0 ± 20.0 


10.0 ± 11.8 


6.8 ± 10.6 


6.2 ± 14.3 


75.1 


3.8 




488 


3.63 ± 0.06 


16.2 ± 22.6 


0.0 ± 25.3 


0.0 ± 38.2 


0.0 ± 34.9 


34.4 ± 25.8 


26.8 ± 21.7 


22.6 ± 18.7 


0.385 




DM Tau 


121 


119.1 ± 2.7 


0.0 ± 15.5 


70.6 ± 18.5 


0.0 ± 13.6 


0.0 ± 9.1 


0.0 ± 6.1 


12.6 ± 7.2 


16.9 ± 8.3 


62.2 


4.5 




260 


0.295 ± 0.040 


0.0 ± 2.4 


25.3 ± 2.7 


34.0 ± 3.2 


40.4 ± 3.3 


0.0 ± 1.7 


0.2 ± 1.7 


0.0 ± 1.2 


2.90 





Table 6 — Continued 









small 


small 


large 


large 








Total 








Temp 




Amorphous 


Amorph. 


Amorph. 


Amorph. 


Crystalline 


Crystalline 


Crystalline 


Dust 






Object 


(K) 


^BB 


Pyroxene^ 


Olivine'' 


Pyroxene'^ 


Olivine'* 


Enstatite" 


Forsterite^ 


Silicas 


Mass 


d.o.f . 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 




DN Tau 


127 


100.1 ± 3.6 


77.1 ± 23.7 


16.4 ± 13.5 


0.0 it 16.2 


0.0 it 10.9 


0.0 it 9.1 


5.9 it 8.1 


0.6 it 12.5 


68.4 


2.9 






431 


1.99 ± 0.03 


0.0 ± 20.5 


49.7 ± 28.9 


0.0 it 31.9 


0.0 it 29.2 


20.6 it 17.0 


7.7 it 13.3 


22.0 it 15.8 


0.283 






DO Tau 


139 


565.9 ± 11.4 


79.4 ± 35.5 


0.0 ± 20.1 


0.0 it 25.1 


0.0 it 16.9 


0.0 it 13.4 


12.7 it 12.8 


7.9 it 15.9 


136 


2.3 






545 


3.59 ± 0.07 


0.0 ± 5.3 


1.3 ± 6.3 


12.1 it 8.2 


76.8 it 11.0 


5.6 it 4.3 


4.2 it 4.0 


0.0 it 2.9 


2.03 






DP Tau 


139 


145.5 ± 4.8 


34.2 ± 6.7 


50.1 ± 7.3 


0.0 it 6.6 


0.0 it 4.4 


4.0 ± 3.6 


4.7 it 3.2 


6.9 it 4.9 


236 


2.9 






488 


2.40 ± 0.05 


0.0 ± 6.0 


64.7 ± 11.8 


0.0 it 10.5 


29.9 it 8.7 


0.0 it 5.1 


5.4 it 4.8 


0.0 it 2.9 


1.16 






DQ Tau 


127 


236.4 dz 7.1 


93.1 ± 30.6 


0.0 ± 15.0 


0.0 it 18.0 


0.0 ± 11.9 


0.0 it 10.6 


6.9 ± 9.9 


0.0 it 13.2 


122 


3.1 






431 


4.43 ± 0.07 


8.1 ± 16.9 


42.8 ± 22.8 


0.0 it 27.6 


0.0 it 25.4 


41.1 it 21.6 


7.9 it 12.1 


0.0 it 9.3 


0.642 






DR Tau 


127 


683.4 ± 17.9 


47.6 ± 10.8 


41.4 ± 9.8 


0.0 ± 10.0 


0.0 it 6.7 


0.0 it 5.5 


3.9 it 4.8 


7.0 it 6.8 


594 


2.1 






545 


6.24 ± 0.10 


0.0 ± 7.4 


59.0 ± 12.5 


5.5 it 11.6 


34.5 it 9.9 


0.5 it 6.0 


0.5 it 5.4 


0.0 it 4.0 


2.32 






DS Tau 


139 


31.1 ± 1.2 


0.0 ± 25.1 


21.9 ± 16.6 


40.8 it 20.2 


0.0 it 13.3 


21.1 it 13.0 


10.9 it 10.6 


5.3 it 13.5 


19.1 


5.1 






545 


0.674 ± 0.014 


0.0 ± 3.7 


45.3 ± 5.0 


44.7 it 5.1 


0.0 it 5.2 


4.1 it 2.9 


3.8 it 2.7 


2.1 it 1.9 


0.682 




1 


F04147+2822 


163 


9.70 ± 0.54 


0.0 ± 8.1 


18.9 ± 5.4 


61.1 it 8.5 


0.0 it 4.9 


4.2 it 3.8 


15.8 it 3.8 


0.0 it 4.7 


24.1 


2.2 


00 




602 


0.361 ± 0.009 


3.2 ± 3.5 


36.2 it 4.1 


60.6 ± 5.4 


0.0 it 4.9 


0.0 it 2.6 


0.0 ± 2.4 


0.0 it 1.9 


0.485 




o 

1 


FM Tau 


127 


93.4 ± 3.0 


60.2 ± 11.1 


27.1 ± 7.9 


0.0 it 9.2 


0.0 it 6.2 


5.6 it 5.2 


3.5 it 4.6 


3.6 it 6.2 


107 


2.4 


1 




431 


1.32 ± 0.03 


7.6 ± 4.4 


78.3 ± 9.3 


10.5 ± 6.4 


0.0 it 6.3 


1.6 it 3.4 


2.0 it 3.1 


0.0 it 2.4 


1.19 






FN Tau 


139 


254.4 it 5.3 


74.6 ± 19.7 


0.0 it 12.9 


0.0 it 14.3 


0.0 it 10.2 


0.1 it 8.6 


25.3 it 9.6 


0.0 it 7.0 


105 


5.1 






488 


1.80 ± 0.04 


0.0 ± 5.8 


38.2 it 7.8 


0.0 it 10.0 


30.6 it 8.3 


28.8 it 6.2 


2.4 it 4.6 


0.0 it 3.2 


1.11 






FO Tau 


115 


133.8 ± 6.0 


43.5 ± 6.6 


23.8 it 5.2 


19.2 ± 5.5 


0.0 it 4.2 


8.9 it 3.5 


4.6 ± 3.2 


0.0 it 4.5 


332 


2.7 






488 


1.57 ± 0.03 


42.6 ± 11.8 


0.0 it 9.6 


40.2 it 12.5 


0.0 it 12.9 


8.1 it 7.0 


2.6 ± 6.7 


6.5 it 5.5 


0.537 






FP Tau 


115 


39.4 ± 2.0 


13.2 ± 6.4 


0.0 ± 5.0 


82.2 it 9.4 


0.0 it 3.9 


1.4 it 3.1 


3.2 ± 2.7 


0.0 ± 4.3 


110 


5.3 






488 


0.400 ± 0.008 


90.2 ± 41.6 


0.0 it 18.4 


0.0 it 27.8 


0.0 it 24.7 


0.9 it 13.1 


8.9 it 11.7 


0.0 it 9.2 


0.0743 






FQ Tau 


151 


12.7 ± 0.4 


47.2 ± 119.2 


0.0 it 95.8 


0.0 ± 116.3 


0.0 ± 80.9 


0.0 it 64.0 


51.7 ± 125.2 


1.1 ± 60.0 


1.02 


3.9 






830 


0.0635 ± 0.0016 


10.0 ± 6.0 


0.0 it 7.3 


4.6 it 10.5 


71.4 it 12.3 


11.9 it 5.2 


1.5 it 4.7 


0.6 it 3.7 


0.0345 






FS Tau 


127 


445.3 ± 12.2 


27.2 ± 9.6 


18.9 it 7.5 


45.4 it 9.8 


0.0 it 6.5 


0.0 it 5.2 


3.6 it 4.5 


4.9 it 6.8 


415 


2.8 






488 


3.05 ± 0.06 


0.0 ± 6.1 


44.8 ± 9.0 


8.2 it 9.9 


40.1 it 9.1 


2.7 it 5.1 


4.3 it 4.6 


0.0 it 3.2 


1.45 






FT Tau 


115 


101.0 ± 4.2 


74.8 ± 18.6 


10.5 ± 11.1 


0.0 it 12.6 


0.0 it 8.2 


0.0 it 7.2 


13.3 it 7.0 


1.4 it 10.5 


108 


5.9 






431 


1.29 ± 0.03 


7.5 ± 6.5 


60.6 ± 11.3 


15.4 it 9.4 


0.0 it 9.7 


7.4 it 5.2 


4.6 it 4.7 


4.4 it 3.8 


0.675 






FV Tau 


127 


359.6 ± 13.4 


35.7 ± 7.2 


58.7 it 8.6 


0.0 it 6.8 


0.0 it 4.5 


1.8 it 3.9 


1.1 it 3.5 


2.7 it 4.8 


652 


3.6 






545 


5.24 ± 0.09 


0.0 ± 3.5 


20.0 it 5.0 


0.0 it 6.8 


72.6 it 8.1 


0.0 ± 3.6 


7.5 it 3.4 


0.0 it 1.5 


3.30 






FX Tau 


115 


58.3 ± 5.4 


0.0 ± 6.5 


55.3 it 6.3 


0.0 it 5.2 


29.8 it 4.1 


1.3 it 3.0 


6.8 it 2.7 


6.7 it 3.8 


358 


5.2 






545 


0.788 ± 0.019 


6.4 ± 3.2 


48.1 it 4.6 


30.3 it 4.3 


0.7 it 4.7 


8.9 it 2.5 


5.5 it 2.4 


0.0 it 1.6 


1.10 







Table 6 — Continued 









small 


small 


large 


large 








Total 






Temp 




Amorphous 


Amorph. 


Amorph. 


Amorph. 


Crystalline 


Crystalline 


Crystalline 


Dust 




Object 


(K) 




Pyroxene^ 


Olivine'' 


Pyroxene*^ 


OUvine'^ 


Enstatite" 


Forsterite^ 


Silicas 


Mass 


2 

X 
d.o.f. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


FZ Tau 


139 


135.1 ± 4.5 


0.0 ± 32.6 


50.9 ± 29.3 


0.0 ± 26.4 


0.0 ± 18.0 


21.6 ± 16.8 


14.7 ± 13.6 


12.8 ± 17.5 


55.7 


7.5 




545 


2.68 ± 0.05 


0.0 ± 6.0 


0.0 ± 6.5 


51.2 ± 9.0 


27.2 ± 7.9 


6.7 ± 4.7 


4.9 ± 4.3 


10.0 ± 3.9 


1.29 




GG Tau A 


139 


192.2 ± 5.0 


73.2 ± 16.0 


8.6 ± 9.7 


0.0 ± 11.9 


0.0 ± 8.1 


0.0 ± 6.5 


12.8 ± 6.6 


5.4 ± 6.9 


130 


2.4 




545 


1.65 ± 0.04 


0.0 ± 2.9 


35.7 ± 3.5 


49.3 ± 4.2 


12.5 ± 3.9 


0.3 ± 2.3 


0.9 ± 2.2 


1.3 ± 1.5 


2.44 




GG Tau B 


139 


46.1 ± 1.0 


51.6 ± 8.7 


37.1 ± 7.2 


0.0 ± 8.0 


0.0 ± 5.5 


0.0 ± 4.3 


8.7 ± 4.2 


2.7 ± 4.8 


38.2 


7.5 




1172 


0.0224 ± 0.0009 


0.0 ± 3.4 


12.6 ± 4.2 


25.5 ± 5.3 


43.1 ± 5.6 


12.7 ± 3.1 


6.2 ± 2.9 


0.0 ± 1.5 


0.0438 




GH Tau 


139 


65.8 ± 2.0 


80.8 ± 15.9 


0.0 ± 9.1 


0.0 ± 10.7 


0.0 ± 7.2 


0.0 ± 5.9 


14.0 ± 5.8 


5.2 ± 7.5 


56.5 


4.3 




716 


0.349 ± 0.007 


0.0 ± 6.3 


1.3 ± 7.3 


10.5 ± 9.6 


76.8 ± 12.8 


7.7 ± 5.0 


0.6 ± 4.6 


3.2 ± 3.5 


0.182 




GI Tau 


115 


171.6 ± 10.1 


0.0 ± 7.6 


0.0 ± 5.7 


87.1 ± 11.0 


0.0 ± 4.0 


9.2 ± 3.5 


3.3 ± 3.2 


0.4 ± 5.4 


610 


3.9 




488 


2.99 ± 0.06 


19.9 ± 3.5 


38.5 ± 4.3 


9.2 ± 5.1 


29.3 ± 4.3 


1.5 ± 2.7 


1.6 ± 2.4 


0.0 ± 1.8 


3.22 




GK Tau 


175 


115.6 ± 2.6 


0.0 ± 10.0 


0.0 ± 7.2 


38.7 ± 8.6 


35.1 ± 7.6 


10.5 ± 4.9 


9.8 ± 4.5 


6.0 ± 5.8 


90.4 


4.2 




1001 


0.283 ± 0.013 


0.0 ± 2.3 


40.3 ± 3.3 


45.6 ± 3.6 


6.6 ± 3.9 


3.8 ± 2.0 


3.4 ± 2.0 


0.3 ± 1.0 


0.992 




GM Aur 


91 


810.0 ± 32.8 


58.4 ± 3.8 


37.7 ± 3.0 


0.0 ± 3.1 


0.0 ± 2.0 


0.0 ± 1.7 


2.1 ± 1.7 


1.9 ± 2.5 


3610 


4.9 




488 


0.326 ± 0.014 


1.0 ± 3.5 


98.1 ± 8.8 


0.0 ± 5.3 


0.0 ± 4.9 


0.0 ± 2.5 


0.9 ± 2.6 


0.0 ± 1.6 


0.714 




ON Tau 


139 


48.9 ± 2.2 


0.0 ± 18.2 


36.4 ± 14.0 


0.0 ± 14.3 


0.0 ± 9.8 


33.6 ± 12.3 


19.9 ± 9.3 


10.1 ± 11.4 


50.7 


5.6 




545 


1.06 ± 0.02 


0.0 ± 3.1 


12.2 ± 3.1 


53.7 ± 4.5 


22.7 ± 3.8 


1.1 ± 2.4 


2.8 ± 2.1 


7.6 ± 1.8 


1.42 




GO Tau 


115 


59.1 ± 1.9 


86.5 ± 17.5 


0.0 ± 9.9 


0.0 ± 10.9 


0.0 ± 7.0 


0.0 ± 6.1 


11.2 ± 5.8 


2.3 ± 8.1 


53.2 


5.3 




545 


0.213 ± 0.006 


0.0 ± 5.1 


34.2 ± 6.0 


42.1 ± 7.1 


14.3 ± 7.1 


4.3 ± 3.9 


3.1 ± 3.3 


2.0 ± 2.8 


0.184 




Haro 6-28 


163 


14.9 ± 0.6 


0.0 ± 13.5 


0.0 ± 9.6 


0.0 ± 12.7 


82.4 ± 19.9 


0.0 ± 6.7 


6.0 ± 5.7 


11.6 ± 7.4 


15.1 


11.8 




659 


0.197 ± 0.006 


11.3 ± 9.5 


0.0 ± 11.1 


44.2 ± 14.8 


0.0 ± 16.4 


21.4 ± 9.0 


14.5 ± 7.8 


8.6 ± 6.1 


0.0825 




Haro 6-37 


163 


74.0 ± 2.0 


61.0 ± 26.4 


1.2 it 18.9 


0.0 ± 22.5 


0.0 ± 15.5 


23.3 ± 14.4 


11.1 ± 11.1 


3.4 ± 13.5 


26.5 


4.9 




716 


0.939 ± 0.022 


6.9 ± 6.5 


28.3 ± 8.0 


0.0 ± 11.6 


36.9 ± 9.8 


18.8 ± 6.0 


3.9 ± 4.9 


5.0 ± 3.8 


0.464 




HK Tau 


115 


360.7 ± 9.4 


91.6 ± 10.0 


0.0 ± 5.1 


0.0 ± 5.7 


0.0 ± 3.7 


0.0 ± 3.4 


8.3 ± 3.1 


0.1 ± 5.0 


506 


6.2 




488 


0.809 ± 0.020 


0.0 ± 6.2 


14.3 ± 7.1 


13.1 ± 9.9 


48.3 ± 9.9 


20.8 ± 5.9 


3.5 ± 4.8 


0.0 ± 3.5 


0.483 




HN Tau 


139 


171.0 ± 5.1 


66.6 ± 12.6 


0.0 ± 8.8 


0.0 ± 9.7 


0.0 ± 6.6 


0.0 ± 5.5 


20.2 ± 5.9 


13.2 ± 6.6 


165 


4.5 




545 


2.01 ± 0.04 


0.0 ± 3.6 


44.6 ± 5.0 


22.8 ± 5.0 


28.1 ± 4.8 


0.0 ± 2.9 


3.0 ± 2.7 


1.5 ± 1.9 


2.07 




HO Tau 


139 


11.7 ± 0.6 


82.7 ± 30.1 


0.0 ± 15.9 


0.0 ± 19.0 


0.0 ± 12.1 


0.0 ± 10.2 


16.9 ± 11.3 


0.4 ± 17.9 


9.81 


5.6 




488 


0.211 ± 0.006 


0.1 ± 5.8 


53.0 ± 8.6 


32.9 ± 7.6 


0.0 ± 8.2 


8.3 ± 4.4 


3.9 ± 3.8 


1.8 ± 3.1 


0.178 




HP Tau 


139 


386.8 ± 9.8 


88.5 ± 31.8 


0.0 ± 15.5 


0.0 ± 18.5 


0.0 ± 12.2 


0.0 ± 10.5 


11.5 ± 10.3 


0.0 ± 18.4 


178 


2.1 




488 


2.98 ± 0.07 


0.0 ± 3.5 


52.8 ± 5.6 


28.4 ± 4.9 


14.6 ± 5.0 


0.0 ± 2.9 


1.8 ± 2.6 


2.3 ± 1.8 


3.55 




HQ Tau 


187 


77.3 ± 2.1 


23.2 ± 5.9 


25.7 ± 4.8 


46.4 ± 6.2 


0.0 ± 4.3 


0.0 ± 3.4 


3.2 ± 3.0 


1.5 ± 4.1 


108 


4.9 




1400 


0.227 ± 0.008 


0.0 ± 2.9 


46.5 ± 4.5 


24.5 ± 4.3 


25.7 ± 4.4 


1.8 ± 2.5 


1.4 ± 2.5 


0.0 ± 1.4 


0.480 





Table 6 — Continued 









small 


small 


large 


large 










Total 






Temp 




Amorphous 


Amorph. 


Amorph. 


Amorph. 


Crystalline 


Crystalline 


Crystalline 


Dust 




Object 


(K) 




Pyroxene* 


Olivine'' 


Pyroxene*^ 


Olivine'' 


Enstatitc" 


Forsterite' 


Silicas 


Mass 


2 

X 

d.o.f. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


IP Tau 


103 


88.6 ± 6.3 


0.0 ± 5.4 


27.7 it 3.9 


64.4 it 6.0 


0.0 it 2.7 


2.6 it 2.5 


0.8 it 2.1 


4.4 it 


4.1 


546 


4.2 




488 


0.648 ± 0.016 


23.7 ± 3.9 


64.9 it 6.8 


5.9 it 5.8 


0.0 it 5.5 


1.4 it 3.1 


4.2 it 2.7 


0.0 it 


1.9 


0.785 




IQ Tau 


139 


51.8 ± 1.9 


92.6 ± 54.1 


0.0 it 29.2 


0.0 it 33.0 


0.0 it 21.9 


0.0 it 17.1 


7.4 it 17.6 


0.0 it 


18.9 


16.9 


6.3 




602 


0.782 ± 0.016 


0.0 ± 4.3 


37.0 ± 5.4 


10.7 ± 6.5 


43.2 it 6.1 


4.2 it 3.5 


4.8 it 3.1 


0.0 ± 


2.3 


0.588 




IS Tau 


175 


13.7 ± 0.5 


0.0 ± 31.1 


0.0 ± 22.0 


0.0 it 26.8 


0.0 it 19.0 


3.8 it 15.2 


64.3 it 36.2 


31.9 ± 


21.4 


5.39 


3.3 




659 


0.417 ± 0.010 


0.0 ± 3.5 


0.0 it 3.7 


41.8 it 4.8 


44.0 it 4.8 


7.1 it 2.7 


2.8 it 2.5 


4.2 it 


2.0 


0.504 




IT Tau 


127 


53.5 ± 2.3 


91.6 ± 17.3 


0.0 it 9.3 


0.0 ± 10.3 


0.0 it 7.0 


3.3 it 5.7 


4.9 ± 5.1 


0.2 it 


7.3 


72.8 


6.3 




659 


0.649 ± 0.011 


0.0 ± 25.3 


37.1 it 30.8 


0.0 ± 40.4 


0.0 it 38.7 


52.6 di 38.6 


10.4 ± 17.6 


0.0 it 


13.8 


0.0709 




Lk Ca 15 


103 


267.6 ± 8.7 


86.5 ± 8.0 


0.0 it 4.5 


0.0 it 4.9 


0.0 it 2.9 


0.0 ± 2.6 


6.2 it 2.8 


7.4 it 


3.8 


562 


7.5 




716 


0.182 ± 0.006 


0.0 ± 2.5 


88.7 it 6.2 


9.7 it 3.8 


0.0 it 3.9 


0.0 it 2.0 


1.6 it 2.0 


0.0 it 


1.1 


0.451 




RW Aur A 


151 


230.2 ± 5.5 


67.2 ± 25.0 


0.0 it 16.0 


0.0 it 19.5 


0.0 it 13.0 


0.0 it 10.8 


18.5 it 11.1 


14.2 it 


14.4 


88.3 


3.9 




602 


2.49 ± 0.05 


4.1 ± 6.2 


68.1 it 12.2 


10.9 it 9.4 


11.9 it 9.0 


3.8 it 5.0 


1.2 it 4.6 


0.0 it 


3.5 


1.29 




UY Aur 


151 


912.9 ± 16.9 


10.3 ± 19.5 


0.0 it 15.2 


72.4 it 25.7 


0.0 it 13.1 


0.0 it 10.1 


4.6 it 8.9 


12.7 it 


11.9 


284 


1.4 




545 


4.96 ± 0.13 


0.0 it 3.6 


59.4 it 6.3 


25.9 it 5.1 


14.7 it 5.3 


0.0 it 2.9 


0.0 it 2.6 


0.0 it 


1.7 


6.23 




UZ Tau E 


175 


93.0 ± 2.4 


47.8 it 10.4 


16.3 ± 7.5 


30.6 it 9.0 


0.0 ± 6.8 


1.0 it 5.3 


1.6 ± 4.6 


2.8 it 


6.8 


77.3 


5.4 




1400 


0.282 ± 0.009 


0.0 it 3.1 


7.0 ± 3.7 


4.8 it 5.3 


85.9 it 7.3 


0.1 it 2.7 


2.3 ± 2.6 


0.0 it 


1.6 


0.401 




V410 Anon 13 


103 


7.28 ± 0.65 


54.4 it 6.8 


6.8 it 5.0 


33.0 it 5.4 


0.0 it 3.6 


1.0 it 3.3 


4.3 it 3.0 


0.4 it 


5.1 


42.5 


5.6 




488 


0.101 ± 0.002 


46.3 it 13.6 


0.0 it 10.5 


10.4 it 14.4 


0.0 it 14.0 


24.0 it 8.9 


13.2 it 7.0 


6.1 it 


6.4 


0.0318 




V710 Tau 


139 


35.5 ± 1.0 


62.9 it 23.1 


0.0 it 15.7 


0.0 ± 18.1 


0.0 ± 11.6 


0.0 ± 10.2 


14.7 ± 10.2 


22.4 ± 


14.8 


16.0 


8.7 




602 


0.427 it 0.008 


0.0 it 8.0 


63.4 it 13.8 


0.0 ± 12.4 


14.1 it 10.3 


14.5 it 6.2 


2.4 it 5.5 


5.6 it 


5.3 


0.176 




V807 Tau 


127 


117.4 ± 3.2 


54.9 ± 7.9 


30.5 ± 6.0 


0.0 it 7.0 


0.0 ± 4.7 


0.0 it 3.7 


8.6 it 3.6 


6.0 it 


4.5 


137 


3.4 




773 


0.428 ± 0.008 


0.0 it 14.6 


8.9 it 16.2 


83.4 it 32.5 


0.0 it 25.1 


2.4 it 12.3 


5.2 it 10.6 


0.0 it 


8.6 


0.0752 




V836 Tau 


127 


37.1 ± 1.4 


58.4 it 16.0 


0.0 it 12.1 


0.0 ± 13.7 


0.0 ± 9.2 


15.8 it 7.9 


16.5 it 7.6 


9.3 ± 


8.6 


32.6 


4.5 




602 


0.237 ± 0.007 


0.0 it 3.3 


24.0 it 3.4 


55.8 it 4.8 


3.7 ± 4.4 


7.8 it 2.4 


4.3 it 2.2 


4.5 it 


1.8 


0.387 




V955 Tau 


115 


81.8 ± 5.9 


30.0 it 8.4 


49.4 it 9.0 


0.0 it 7.8 


0.0 it 5.2 


3.5 it 4.6 


15.4 it 4.7 


1.7 it 


6.5 


257 


2.5 




488 


2.06 ± 0.03 


12.3 it 6.1 


0.0 it 6.5 


63.0 it 10.0 


0.0 it 8.6 


11.5 it 4.8 


4.8 ± 4.3 


8.4 ± 


3.9 


0.966 




VY Tau 


115 


47.7 ± 1.7 


62.9 it 38.7 


0.0 it 30.0 


0.0 ± 31.3 


0.0 ± 22.0 


0.0 ± 17.1 


37.1 ± 27.4 


0.0 it 


15.5 


15.3 


7.2 


LLtl 


545 


0.226 ± 0.006 


22.0 it 6.2 


19.7 it 6.5 


43.8 it 8.7 


0.0 ± 9.1 


5.0 ± 4.8 


8.0 it 4.5 


1.5 it 


3.6 


0.153 




XZ Tau 


151 


787.4 ± 18.1 


75.3 it 22.1 


20.5 it 12.5 


0.0 it 16.8 


0.0 it 12.2 


0.0 it 8.3 


4.2 it 7.8 


0.0 it 


8.9 


317 


2.0 




488 


12.1 ± 0.2 


0.0 it 9.6 


13.4 it 10.4 


68.3 it 17.7 


10.4 it 14.4 


4.8 ± 7.7 


3.0 ± 7.0 


0.0 it 


5.5 


3.54 




TIL Tau 


211 


6.43 ± 0.15 


0.0 it 33.3 


49.4 it 34.8 


0.0 ± 35.2 


0.0 ± 28.2 


0.4 ± 18.9 


0.0 ± 15.6 


50.1 it 


35.1 


1.17 


3.5 




1400 


0.0106 ± 0.0007 


8.6 it 4.3 


0.0 it 5.4 


78.9 it 10.1 


0.0 it 8.5 


7.0 it 3.8 


1.2 ± 3.7 


4.3 it 


2.4 


0.0218 
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small 


small 


large 


large 








Total 






Temp 




Amorphous 


Amorph. 


Amorph. 


Amorph. 


Crystalline 


Crystalline 


Crystalline 


Dust 




Object 


(K) 




Pyroxene^ 


Olivine'' 


Pyroxene'^ 


Olivine'^ 


Enstatite" 


Forsterite^ 


Silicas 


Mass 


d.o.f. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


ZZ Tau IRS 


127 


280.3 ± 7.5 


84.4 ± 11.1 


11.0 ± 5.9 


0.0 ± 7.3 


0.0 ± 4.8 


0.0 ± 3.7 


4.6 ± 3.9 


0.0 ± 4.7 


323 


4.9 




545 


1.28 ± 0.03 


0.0 ± 3.1 


0.0 ± 4.1 


0.0 ± 6.2 


96.3 ± 8.8 


0.0 ± 2.9 


3.7 ± 2.8 


0.0 ± 1.5 


1.18 





Note. — Col. (1): Object name. Col. (2): One of two dust model temperatures (Kelvin) . Col. (3): Solid angle, fiss, of blackbody of temperature specified 
in Col. (2) representing continuum emission, expressed in units of 10^^^ steradians. Cols. (4)-(10): mass percentages and their propagated uncertainties 
of all dust masses at temperature specified in Col. (2). One dust model is completely specified by two adjacent rows - the row following the object's name 
and the row beneath that one. Propagated uncertainties are the square root of the sum of the squares of terms, such that each term is the product of the 
uncertainty in one of seven mass weights and the partial derivative of the mass fraction of the dust species in question with respect to the same rr^s weight. 
This way, even if a species has zero mass, its uncertainty still contributes to the propagated uncertainty in question. Col. (11): Total dust riiSss at one 
temperature in 10~* lunar masses, computed assuming 140pc to Taurus-Auriga. Col. (12): P^r degree of freedom, determined over 7.7 < A < jil/im. 



'^Optical constants for amorphous pyroxene of cosmic composition from I Jaeger et al.l (|l994l ). assuming CDE2 (jFabian et al.ll2001l ) 
''Optical constants for amorphous olivine MgFeSi04 from lDorschner et al.l (jl995[ ). assuming CDE2 

'^Optical constants for amorphous pyroxene of cosmic composition from I Jaeger et al" (1994), using the Bruggeman EMT and Mie theory (Bohren & Huffman 
1983) with a volume fraction of vacuum oi f — 0.6 for porous spherical grains of radius 5 /im 

''Optical constants for amorphous olivine MgFeSi04 from lDorschner et al.l (jl995l ). using the Bruggeman EMT and Mie theory (Bohren & Huffman 1983) 
with a volume fraction of vacuum of / = 0.6 for porous spherical grains of radius 5 /Ltm 

"Opacities for clinoenstatite Mgg.gFeo.iSiOs from lchihara et aL (2002) 



^Optical constants for 3 crystallographic axes of forsterite, Mg2Si04, from lSogawa et al.l (j2006l ). assuming tCDE (see discussion in Section 4.2.3). 
s Opacity for annealed sihca by Fabian et al. ( 200Cll ). 
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Table Al: Mispointing Corrections 

LL2nodl or LL2nod.2 or LLlnodl or LLlnod2 or 



Object 


modules 


SL2nodl 


SL2nod2 


SLlnodl 


SLlnod2 


SHnodl 


SHnod2 


LHnodl 


LHnoc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


AA Tau 


SLLL 


1.00 


1.00 


1.00 


1.00 


1.02 


1.00 


1.02 


1.04 


BP Tau 


SLLL 


1.01 


1.00 


1.05 


1.00 


1.00 


1.01 


1.00 


1.01 


CI Tau 


SLLL 


1.00 


1.00 


1.00 


1.01 


1.00 


1.03 


1.00 


1.01 


CoKu Tau/3 


SLLL 


1.01 


1.00 


1.00 


1.00 


0.96 


1.00 


1.00 


1.00 


CoKu Tau/4 


SLLL 


1.01 


1.00 


1.00 


1.01 


1.00 


1.00 


1.00 


1.01 


CW Tau 


SLSHLH 


1.13 


1.14 


1.10 


1.09 


1.25 


1.00 


1.00 


1.00 


CX Tau 


SLLL 


1.12 


1.00 


1.24 


1.20 


1.21 


1.20 


1.20 


1.21 


CY Tau 


SLLL 


1.01 


1.00 


1.00 


1.01 


1.01 


1.00 


1.00 


1.00 


DD Tau 


SLLL 


1.00 


1.00 


1.01 


1.00 


1.00 


1.01 


1.00 


1.00 


DE Tau 


SLLL 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.01 


DP Tau 


SLSHLH 


1.00 


1.00 


1.00 


1.00 


1.17 


0.96 


1.02 


1.00 


DG Tau 


SLSHLH 


1.00 


1.00 


1.05 


1.00 


1.12 


1.00 


1.05 


1.06 


DH Tau 


SLLL 


1.00 


1.00 


1.00 


1.00 


1.00 


1.01 


1.00 


1.03 


DK Tau 


SLSHLH 


1.01 


1.00 


1.00 


1.01 


1.00 


1.01 


1.00 


1.00 


DL Tau 


SLSHLH 


1.01 


1.00 


1.00 


1.00 


1.00 


1.04 


1.00 


1.01 


DM Tau 


SLLL 


1.15 


1.15 


1.15 


1.15 


1.00 


1.00 


1.00 


1.01 


DN Tau 


SLLL 


1.00 


1.00 


1.00 


1.00 


1.02 


1.00 


1.00 


1.01 


DO Tau 


SLSHLH 


1.01 


1.00 


1.00 


1.00 


1.00 


1.07 


1.05 


1.06 


DP Tau 


SLLL 


1.00 


1.00 


1.01 


1.00 


1.01 


1.00 


1.00 


1.01 


DQ Tau 


SLLL 


1.05 


1.05 


1.05 


1.05 


1.00 


1.00 


1.00 


1.00 


DR Tau 


SLSHLH 


1.00 


1.00 


1.00 


1.00 


1.05 


0.95 


1.01 


1.00 


DS Tau 


SLLL 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


P04147+2822 


SLLL 


1.00 


1.00 


1.00 


1.01 


1.00 


1.01 


1.01 


1.00 


PM Tau 


SLLL 


1.02 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


PN Tau 


SLLL 


1.02 


1.00 


1.10 


1.00 


1.10 


1.10 


1.10 


1.10 


PO Tau 


SLLL 


1.03 


1.00 


1.15 


0.92 


1.00 


1.01 


1.00 


1.00 


PP Tau 


SLLL 


1.01 


1.00 


1.00 


1.01 


1.00 


1.00 


1.00 


1.00 


PQ Tau 


SLLL 


1.00 


1.01 


1.00 


1.00 


1.01 


1.00 


1.00 


1.00 


PS Tau 


SLSHLH 


1.00 


1.00 


1.02 


1.00 


1.09 


0.93 


1.00 


1.03 


PT Tau 


SLLL 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


PV Tau 


SLSHLH 


1.01 


1.00 


1.04 


1.00 


1.12 


0.97 


1.00 


1.00 


PX Tau 


SLLL 


1.01 


1.00 


1.02 


1.00 


1.00 


1.00 


1.05 


1.06 


PZ Tau 


SLSHLH 


1.01 


1.00 


1.03 


1.00 


1.19 


0.97 


1.01 


1.00 


GG Tau A 


SLLL 


1.04 


1.03 


1.04 


1.03 


1.00 


1.00 


1.00 


1.00 


GG Tau B 


SLLL 


1.10 


1.15 


1.10 


1.15 


1.02 


1.00 


1.01 


1.00 


GH Tau 


SLLL 


1.00 


1.00 


1.00 


1.01 


1.00 


1.01 


1.00 


1.01 


GI Tau 


SLLL 


1.02 


1.00 


1.04 


1.00 


1.00 


1.00 


1.00 


1.01 


GK Tau 


SLSHLH 


1.02 


1.00 


1.03 


1.00 


1.06 


1.05 


1.05 


1.06 


GM Aur 


SLLL 


1.00 


1.00 


1.00 


1.02 


1.00 


1.01 


1.00 


1.00 


GN Tau 


SLLL 


1.01 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


GO Tau 


SLLL 


1.03 


1.00 


1.00 


1.00 


1.00 


1.01 


1.00 


1.01 


Haro 6-28 


SLSHLH 


1.03 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


Haro 6-37 


SLSHLH 


1.00 


1.00 


1.03 


1.00 


1.06 


1.00 


1.00 


1.00 


HK Tau 


SLLL 


1.00 


1.00 


1.00 


1.01 


1.00 


1.01 


1.00 


1.01 


HN Tau 


SLSHLH 


1.01 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


HO Tau 


SLLL 


1.01 


1.00 


1.00 


1.00 


1.06 


1.00 


1.02 


1.00 


HP Tau 


SLSHLH 


1.02 


1.00 


1.00 


1.01 


1.07 


1.13 


1.06 


1.10 


HQ Tau 


SLSHLH 


1.01 


1.00 


1.00 


1.01 


1.00 


1.11 


1.03 


1.03 


IP Tau 


SLLL 


1.01 


1.00 


1.01 


1.00 


1.00 


1.00 


1.00 


1.01 


IQ Tau 


SLLL 


1.01 


1.00 


1.00 


1.00 


1.02 


1.04 


1.02 


1.04 


IS Tau 


SLLL 


1.00 


1.02 


1.00 


1.02 


1.00 


1.00 


1.00 


1.00 


IT Tau 


SLLL 


1.02 


1.00 


1.00 


1.04 


1.00 


1.01 


1.00 


1.01 


LkCa 15 


SLLL 


1.01 


1.00 


1.00 


1.02 


1.00 


1.02 


1.00 


1.00 


RW Aur A 


SLSHLH 


1.00 


1.00 


1.00 


1.00 


1.01 


1.00 


1.01 


1.00 


UY Aur 


SLSHLH 


1.00 


1.00 


1.00 


1.00 


1.11 


1.01 


1.08 


1.07 


UZ Tau E 


SLSHLH 


1.03 


1.00 


1.00 


1.02 


1.37 


0.87 


1.10 


1.11 


V410 Anon 13 


SLLL 


1.02 


1.00 


1.00 


1.02 


1.08 


1.00 


1.05 


1.00 


V710 Tau 


SLLL 


1.02 


1.00 


1.12 


1.07 


1.04 


1.00 


1.05 


1.06 


V807 Tau 


SLLL 


1.01 


1.00 


1.00 


1.01 


1.01 


1.00 


1.00 


1.00 


V836 Tau 


SLLL 


1.00 


1.00 


1.00 


1.01 


1.00 


1.00 


1.00 


1.00 


V955 Tau 


SLLL 


1.00 


1.00 


1.00 


1.01 


1.00 


1.00 


1.00 


1.01 


VY Tau 


SLLL 


1.00 


1.00 


1.01 


1.00 


1.01 


1.00 


1.00 


1.01 


XZ Tau 


SLSHLH 


1.00 


1.01 


1.01 


1.00 


1.00 


1.01 


1.00 


1.02 


ZZ Tau" 


SLLL 


















ZZ Tau IRS 


SLLL 


1.00 


1.00 


1.03 


1.00 


1.00 


1.00 


1.00 


1.00 



Note. — Col. (1): Object name. Cols. (2)-(9): Multiplicative scalars applied to the spectrum of one 
order of one nod for a given object to match the flux density of the other nod as described in the text. SL2 
means Short-Low order 2, SLl is Short-Low order 1, LL2 is Long-Low order 2, LLl is Long-Low order 1. 
Tor ZZ Tau, there were 6 DCEs per spectral order. Each DCE of each order was normalized to the DCE 
with the highest average flux density over that order, except for Long-Low Bonus order. In that case, the 
DCEs were n ormalized to the av erage flux level of the DCE with the second-highest flux density level (see 
discussion by Sargent et al. 2008h . 



